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ABSTRACT 


This  thesis  devoted  to  studying  the  properties  of  a  digital  M-ary  signaling  scheme  called 
MFSK'QPSK  which  combines  the  features  of  M-ary  FSK  and  QPSK  modulation.  This 
scheme  is  shown  to  be  a  form  of  biorthogonal  signaling,  and  a  receiver  that  is  optimum 
in  error  rate  sense  is  introduced  and  analyzed  in  order  to  obtain  its  bit  error  rate  (BER) 
performance  when  signal  interference  can  be  modeled  as  additive  white  Gaussian  noise. 
The  bandwidth  of  the  MFSK/QPSK  signaling  scheme  is  derived  and  compared  with  that 
of  M-ary  PSK  and  M-ary  FSK  signaling.  MFSK/QPSK  signaling  results  in  receivers 
that  provide  error  rates  much  lower  than  corresponding  receivers  for  M-ary  PSK  and 
slightly  lower  than  those  for  M-ary  FSK,  while  maintaining  a  bandwidth  efficiency  at 
least  three  times  greater  than  that  of  M-ary  FSK  though  not  as  great  as  that  of  M-ary 
PSK  signaling.  A  significant  effort  has  been  devoted  to  examining  a  detection  method 
known  as  direct  bit  detection  (DBD)  which  has  been  applied  to  M-ary  PSK  and 
MFSK  QPSK  signaling.  DBD  receivers  are  introduced,  and  the  BER  of  the 
MFSK  QPSK  DBD  receiver  is  derived  and  compared  with  that  of  the  optimum  (symbol 
discrimination)  receiver.  Lastly,  the  effects  on  performance  due  to  Rayleigh  fading  are 
analyzed  and  evaluated  for  M-ary  PSK  DBD  receivers  as  well  as  for  the  MFSK'QPSK 
optimum  (symbol  discrimination)  and  DBD  receiver. 
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I.  INTRODUCTION 


A.  BACKGROUND 

There  exists  an  increasing  need  for  high-data-rate,  reliable  digital  communication 
systems.  The  reasons  for  this  need  stem  from  the  desire  by  an  increased  number  of  users 
to  transmit  larger  amounts  of  data  via  digital  communication  schemes  which  provide 
more  reliability,  flexibility,  and  lower  cost  w’hen  compared  with  analog  communication 
systems.  These  increased  demands  come  at  a  time  when  the  available  bandwidth  in  the 
free  space  channel  is  becoming  a  limited  resource.  These  data  transmission  requirements 
coupled  with  limited  channel  resources  are  the  driving  force  behind  the  continual  search 
for  digital  modulation  techniques  which  operate  at  low  error  rates  and  are  bandwidth 
efficient.  Present  signaling  schemes  in  use  are  the  well  known  M-ary  Phase  Shift  Keying 
(M-PSK)  and  M-ary  Frequency  Shift  Keying  (M-FSK).  These  signaling  schemes  are 
briefly  described  in  Sections  C  and  D  of  this  chapter,  and  the  corresponding  receiver 
structures  for  recovering  the  digital  information  are  described  in  Section  E,  while  the 
performance  of  these  receivers,  specified  in  terms  of  the  symbol  error  rates,  for  both 
M-ary  PSK  and  M-ary  FSK  modulation  are  given  in  Section  F.  The  bandwidth  re¬ 
quirements  associated  with  these  two  modulation  schemes  are  given  in  Section  G. 

B.  CONTENTS  AND  SUMMARY  OF  RESULTS 

This  thesis  explores  a  hybrid  form  of  digital  modulation,  namely  MFSK  QPSK, 
which  embodies  features  of  both  M-ary  FSK  and  Quaternary  (or  four  phase)  PSK.  The 
analysis  compares  and  contrasts  the  properties  of  this  modulation  scheme  to  the  widely 
known  M-ary  PSK  and  M-ary  FSK  modulation.  In  Chapter  3,  MFSK, QPSK  is  shown 
to  be  a  form  of  biorthogonal  signaling,  and  therefore  its  corresponding  receiver  per¬ 
formance  is  easily  revealed  given  what  is  already  known  for  this  class  of  signals.  The 
optimum  MFSK'QPSK  receiver  that  recovers  the  digital  information  with  minimum 
probability  of  error  under  specific  noise  interference  assumptions  is  analyzed  from  a 
performance  standpoint,  which  traditionally  has  been  given  in  terms  of  the  Symbol  Error 
Rate  (SER).  Because  a  more  meaningful  performance  parameter  is  the  receiver's  Bit 
Error  Rate  (BER),  Chapter  3  has  been  devoted  to  deriving  the  receiver's  BER  as  a 
function  of  signal-to-noise  ratio,  along  with  the  bandwidth  requirements  and  the  band¬ 
width  efficiency  of  the  MFSK, QPSK  scheme.  Comparisons  are  made  between  the 
MFSK  QPSK  scheme  and  conventional  phase  and  frequency  modulation  methods  with 
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respect  to  SER  performance  and  bandwidth  efficiency  in  Chapter  5.  The  results  show 
that  the  MFSK, QPSK  scheme  when  compared  to  M*ary  PSK,  does  not  attain  the  high 
bandwidth  efficiencies  achieved  by  M-ary  PSK,  yet  outperforms  M-ary  PSK  substan¬ 
tially  in  terms  of  SER  performance,  while  performing  slightly  better  than  M-ary  FSK  in 
terms  of  SER  performance  and  only  requiring  one  third  to  one  fourth  the  bandwidth 
required  by  an  M-ary  FSK  system. 

Additionally,  this  thesis  addresses  the  properties  of  a  class  of  digital  communication 
receivers  which  decode  the  individual  bits  that  make  up  a  symbol  directly,  without  need 
for  symbol-to-bit  mapping  post-detection  circuitry.  These  are  referred  to  as  direct  bit 
detection  (DBD)  receivers,  because  each  of  the  bits  that  make  up  a  symbol  are  deter¬ 
mined  independently.  These  receivers  have  been  proposed  and  analyzed  for  M-ary  PSK 
Ref.  1  [pp.  461-463)  as  well  as  for  the  MFSK  QPSK  signaling  scheme  Ref.  2  [pp.  16-101] 
under  discussion  here.  DBD  receiver  structures  for  M-ary  PSK  are  introduced  in 
Chapter  2,  and  previously  derived  expressions  for  their  BER  performance  are  presented. 
In  Chapter  3,  a  DBD  receiver  structure  for  MFSK/QPSK  signaling  is  described,  and 
analysis  is  carried  out  in  order  to  derive  its  BER  performance,  with  corresponding  results 
presented  in  graphical  form  in  Chapter  5.  Also  in  Chapter  5  the  BER  performance  of 
the  MFSK  QPSK  DBD  receiver  is  compared  to  that  of  the  optimum  receiver  that  per¬ 
forms  symbol  detection  followed  by  symbol-to-bit  mappings  prior  to  delivering  binary 
data  to  the  intended  user.  Furthermore,  the  results  demonstrate  that  while  the  DBD 
receiver  does  not  perform  as  well  as  the  optimum  receiver  in  BER  sense,  its  simplified 
receiver  structure  and  direct  bit  detection  properties  offer  other  potential  advantages. 

In  choosing  a  modulation  method,  it  is  important  that  the  effect  of  system  disturb¬ 
ances  such  as  Rayleigh  fading,  which  is  common  in  digital  microwave  radio  transmission 
applications,  be  analyzed  in  order  to  fully  understand  all  aspects  of  expected  system 
performance  for  the  different  signaling  schemes  under  consideration.  The  effect  of 
Rayleigh  fading  on  the  BER  performance  of  the  MFSK  QPSK  receivers  as  well  as  of  the 
M-ary  PSK  DBD  receivers  is  analyzed  and  evaluated  in  Chapter  4.  The  results  demon¬ 
strate  that  in  all  cases,  Rayleigh  fading  can  severely  degrade  receiver  performance. 

C.  M-ARY  PHASE-SHIFT  KEYING 

M-ary  signaling  techniques  involve  forming  groups  of  k  bits  to  generate  M  =  2* 
symbols  which  are  assigned  to  M  different  waveforms  suitable  for  transmission  over  a 
channel.  One  form  of  M-ary  signaling,  namely  M-ary  Phase-Shift  Keying  (M-PSK),  is 
a  bandwidth  efficient  digital  modulation  scheme.  The  phase  of  the  transmitted 
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sinusoidal  carrier  waveform  can  take  on  one  of  M  different  phases  to  represent  M  dif¬ 
ferent  symbols.  The  M-arv  PSK  waveform  can  be  described  mathematically  as 


¥')  = 


(1.1) 


for  i  =  0,  1,...,  M-l  and  0  <,  t^T„  where  generally  (as  is  the  case  here)  the  phase  dif¬ 
ference  between  waveform  pairs  is  taken  to  be  constant.  Each  symbol  or  waveform 
sit)  is  transmitted  at  a  rate  of  R,  =  -y-  symbols/ sec,  where  T,  is  the  time  duration  of  each 
waveform  s,(r).  For  k  =  2  and  M  =  4,  the  resulting  modulation  form  is  known  as 
Quaternary  Phase-Shift  Keying  (QPSK),  where  the  phase  difference  between  waveform 
pairs  is  ■—  radians.  For  8-PSK,  in  which  k  —  3  and  M  =  8,  waveform  pairs  are  separated 
by  a  phase  of  -y  radians.  Signal  constellations  for  QPSK  and  8-PSK  as  well  as  symbol 
assignments  for  each  signal  (or  phasor  representation)  are  presented  in  Figure  1  on  page 
4. 


D.  M-ARY  FREQUENCY-SHIFT  KEYING 


M-ary  Frequency- Shift  Keying  (M-FSK)  is  another  form  of  M-ary  signaling,  in 
which  the  waveforms  are  mathematically  expressed  as 


cos  2 nft 


(1.2) 


where  the  frequencies/,  typically  are  given  by 

fi  =/>  +  'A/ 


(1-3) 


for  i  =  0,  1,....  M-l  and  0  <,  t  <,Tr  The  signal  duration  is  normally  chosen  so  that  fT, 
is  an  integer  or  at  least  fT,  >  >  1  holds.  The  frequency  separation  A/ is  often  set  equal 
to  R,  providing  minimum  frequency  spacing  between  signals  s,(/)  ,  while  forming  a  set 
of  orthogonal  signals.  For  k  =  2  and  M  =  4,  the  resulting  modulation  is  Quadrature 
Frequency-Shift  Keying  (QFSK).  Magnitude  spectra  for  QFSK  and  M-ary  FSK  (  for 
arbitrary'  M  )  are  shown  in  Figure  2  on  page  5  and  in  Figure  3  on  page  6,  respectively. 

E.  RECEIVERS  OF  M-ARY  PSK  AND  M-ARY  FSK 

Signals  transmitted  through  a  physical  channel  always  encounter  interference.  The 
SER  results  given  later  in  this  chapter  and  the  BER  analysis  presented  in  Chapter  3  are 
obtained  under  the  assumption  that  this  interference  can  be  modeled  as  a  sample  func¬ 
tion  of  an  additive  white  Gaussian  noise  (AWGN)  process  having  double  sided  power 
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spectral  density  (PSD)  of-y-.  L'nder  this  assumption,  receivers  can  be  derived  that  are 
optimum  in  the  sense  that  the  probability  of  a  symbol  error  is  minimized.  One  structural 
form  of  an  optimum  receiver  for  M-ary  PSK  is  shown  in  Figure  4  on  page  8.  This  re¬ 
ceiver  as  well  as  all  receivers  considered  herein  is  assumed  to  be  coherent,  meaning  that 
the  phase  of  the  sinusoidal  carrier  is  known  at  the  receiver.  The  concept  of  coherence 
applies  to  this  receiver  in  that  the  cosine  and  sine  functions  which  are  correlated  with  the 
received  signal  are  in  phase  with  the  unmodulated  carrier  sinusoid.  Because  the  receiver 
correlates  the  received  signal  with  a  set  of  sinusoidal  functions,  which  in  this  case  con¬ 
sists  of  each  of  the  possible  transmitted  signals,  it  is  referred  to  as  a  correlation  receiver. 

An  optimum  receiver  structure  for  M-ary  FSK  is  shown  in  Figure  5  on  page  9, 
which  also  corresponds  to  a  correlation  receiver.  The  correlator  with  a  local  signal  input 
which  most  closely  matches  the  incoming  signal,  will  produce  in  the  absence  of  interfer¬ 
ence,  the  largest  output  so  that  the  receiver  can  decide  which  signal  was  sent  according 
to  which  correlator  output  is  the  largest. 

F.  SYMBOL  ERROR  RATES  FOR  M-ARY  PSK  AND  M-ARY  FSK 

For  coherent  demodulation  of  M-ary  PSK  using  the  optimum  receiver  shown  in 
Figure  4,  the  SER  is  given  by  Ref.  3  (pp.  204-207] 
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where 


£(£)  = 


(1-5) 


The  integral  on  the  right  hand  side  of  Eq.  1.4  has  no  known  closed  form,  and  therefore 

£ 

must  be  evaluated  numerically  if  plotting  of  P,{e)  as  a  function  of  -rr-  is  desired.  For 

E  ’<> 

high  values  of  — ,  the  SER  can  be  tightly  upperbounded  by 


PAe) 


<20 


(1.6) 
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Figure  4. 


W> 


Receiver  Structure  for  M-ary  PSK 


8 


Figure  5. 


cos(2  n 


Receiver  Structure  for  M-ary  FSK 


9 


For  M-ary  FSK  modulation  and  an  AWGN  interference  model,  the  optimum  co¬ 
herent  receiver  of  Figure  5  on  page  9  achieves  a  probability  of  correctly  decoding  a 
symbol  given  by  Ref.  3  [pp.  200-203] 


ps{Q  = 

*oo 

f-v+ V-T 

I  g(x)dx 

* 

— oo 

— OO 

-\M-\ 


g(y)dy 


(1.7) 


Ao 


where  a  =  -j-,  and  g(Q  is  the  zero  mean  unit  variance  Gaussian  distribution,  or 

g( 0  =  — F="  expf  -  4" 

\  1 


(1.8) 


This  expression  can  be  rewritten  as 


(1-9) 


£ 

Also,  the  SER  can  be  closely  upperbounded  for  high  values  of  -rr-  by 

/>,{*}  <(M-1 


(1.10) 


G.  BANDWIDTH  REQUIREMENTS 

For  M-ary  PSK,  the  carrier  is  switched  at  the  rate  of  Rs  times  per  second.  Since 
there  are  k  bits  associated  with  each  symbol,  the  bit  rate  Rb  can  be  expressed  as 

Rb  =  kRs  (1.11) 

The  null-to-null  bandwidth  W  for  M-ary  PSK  is  given  by 


IV  =  2RS 


(1.12) 


For  M-ary  FSK  the  carrier  is  also  switched  at  the  symbol  rate  R, .  Looking  at  Fig¬ 
ure  3  it  can  be  seen  that  the  total  null-to-null  bandwidth  left-most  main  lobe  to  right¬ 
most  main  lobe  for  M-ary  FSK  can  be  written  as 
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(1.13) 


W  =  (M  4-  \)RS 


(2k+l)Rb 

k 


From  Eq.  1.13  it  is  apparent  that  the  null-to-null  bandwidth  of  M-ary  FSK  increases 
with  increasing  values  of  k,  however  Eq.  1.12  shows  that  the  null-to-null  bandwidth  of 
M-ary  PSK  decreases  with  increasing  k.  This  property  makes  M-ary  PSK  a  desirable 
modulation  method  when  bandwidth  must  be  minimized.  In  the  next  chapter  M-ary 
PSK  receivers  of  a  different  structure  than  the  one  shown  in  Figure  4  are  introduced, 
and  their  bit  error  rates  presented. 
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II.  DIRECT  BIT  DETECTION  METHODS  FOR  M-ARY  PSK 

A.  M-ARY  PSK  DIRECT  BIT  DETECTION  RECEIVERS 

M-ary  PSK  modulation  as  discussed  in  Chapter  1,  Section  A,  is  a  signaling  scheme 
that  achieves  bandwidth  conservation  at  the  expense  of  signal-to-noise  ratio  required  to 
maintain  a  certain  bit  error  rate  performance  level.  The  two  forms  of  M-ary  PSK  sig¬ 
naling  discussed  in  this  chapter  are  8-PSK  and  16-PSK.  For  8-PSK,  k  =  3  bits/symbol 
and  M  —  8  symbols,  while  for  16-PSK,  k  —  4  bits/symbol  and  M=  16  symbols.  The 
M-ary  PSK  receiver  discussed  in  Chapter  1,  Section  E  and  most  other  optimum  receiver 
structures  discussed  in  the  literature,  require  symbol-to-bit  mapping  to  be  carried  out 
after  a  decision  has  been  made  as  to  what  signal  was  transmitted.  A  receiver  structure 
which  does  not  require  this  symbol-to-bit  mapping  is  proposed  in  Ref.  4  [pp.  3-5]  and  in 
Ref.  1  [pp.  461-463]  for  8-PSK  and  16-PSK,  respectively.  Block  diagrams  for  these  re¬ 
ceivers  are  shown  in  Figure  6  on  page  13  for  8-PSK,  and  in  Figure  7  on  page  14  for 
16-PSK.  Because  each  of  the  bits  are  determined  independently  by  these  receivers  by 
appropriate  processing  of  the  correlator  outputs,  these  receivers  are  referred  to  as  direct 
bit  detection  (DBD)  receivers.  The  processes  used  by  DBD  receivers  for  recovering  bits 
make  it  possible  to  adjust  the  BER  of  the  individual  bits  within  a  symbol  by  means  of 
varying  the  angles  between  the  signal  waveforms.  These  angle  adjustments  are  shown 
in  Figure  8  on  page  15  for  8-PSK,  and  in  Figure  9  on  page  16  for  16-PSK,  where  a  and 
/?  are  the  adjustable  phase  parameters. 

B.  PERFORMANCE  OF  DBD  RECEIVERS  FOR  M-ARY  PSK 

A  BER  analysis  for  the  8-PSK  DBD  receiver  was  carried  out  by  Ref.  4  [pp.  5-14]  in 
which  a  coherent  receiver  was  assumed,  and  AWGN  used  as  the  interference  model. 
The  three  bits  associated  with  each  symbol  are  identified  as  IB,  2B  and  3B  corresponding 
to  the  most  significant  bit,  the  second-most  significant  bit  and  the  least  significant  bit, 
respectively.  The  results  of  the  BER  analysis  showed  that  the  individual  bit  error  prob¬ 
abilities  were  not  equal  for  the  IB,  2B  and  3B.  In  fact  these  probabilities  are  given  by 

JV{1B  incorrect}  =  ■}  cos  a)  +  \  Q sina)  (J  () 

=  Pr{3B  incorrect} 


and 
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Figure  6.  Block  Diagram  of  the  8-PSK  DBD  Receiver 
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Figure  7.  Block  Diagram  of  the  16-PSK  DBD  Receiver 


sin(27tfoO 


Figure  8.  Signal  Constellation  for  8-PSK  Showing  Variable  Phase  Parameters 
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lv 


Showing  Variable  Phase  Parameters 


(2.2) 


(2.3) 


For  16-PSK,  there  are  four  bits  associated  with  each  symbol  which  are  identified  as 
IB,  2B,  3B  and  4B  corresponding  to  the  most  significant  bit,  the  second-most  significant 
bit,  the  third-most  significant  bit,  and  the  least  significant  bit,  respectively.  The  error 
probability  associated  with  each  of  these  bits  is  different  as  given  by  Ref.  1  [pp.  461-463], 
under  the  assumption  of  a  coherent  receiver  and  an  AWGN  interference  model.  The 
bit  error  probabilities  for  the  individual  bits  in  16-PSK  are  given  by 


Pr{  1 B  incorrect)  = 


=  /V{4B  incorrect) 


Pr{ 2B  incorrect)  =  -y  '[  cos  a  -  sin  a]j  4-  Q^yJ  ~y~  [  cos  «  +  sin  a]'j  • 

^"2(V^tcos“_sina]) 

L  _  J  _  (2.‘ 

+  Y  [cos  p-  sin  fi"]j  +  [  cos  p  +  sin  0]J  . 

Y  ~  Q^yf-^l  cos  p  -  sin 


[cos  P  -  sin /?] 
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/V{3B  incorrect}  = 


4 


where  the  integral  in  Eq.  2.6  has  not  been  solved  analytically,  but  has  been  evaluated 
numerically  in  Ref.  1  [pp.  466-467],  and  the  functions  /(I-,  a,  /?)  and  g,(T,  «,  p)  are  de¬ 
fined  in  Table  I. 


Table  1,  FUNCTION  DEFINITIONS 


Index 

i 

f{ r, P) 

i,(r,  a,  P) 

1 

sin(T  —  a) 

cos(T  —  a) 

2 

cos(T  -  a) 

— 

a 

1 

U 

<Sj 

3 

cos(T  4-  a) 

sin(T  4-  a) 

4 

sin(T  +  a) 

cos(T  +  a) 

5 

sin(T  -  /?) 

cos(T  —  p) 

6 

cos(T  —  P) 

sin(r  —  p) 

7 

cos(T  4-  P) 

sin(T  +  P) 

8 

sin(E  +  P) 

cos(f  4-  P) 

The  overall  BER  for  16-PSK  can  then  be  expressed  as 

Pb{e)  = 

4-  /V{1B  incorrect}  +  4-  /V{2B  incorrect}  +  -7  /V{'3 B  incorrect} 

L  A  4 
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These  results  have  been  presented  here  so  that  results  in  the  sequel  can  be  compared  to 
the  BER  results  associated  with  DBD  coherent  receivers  for  8-PSK  and  16-PSK  signal¬ 
ing. 
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III.  MFSK/QPSK 


A.  MFSK/QPSK  AS  A  BIORTHOGONAL  SIGNAL  SET 

If  the  features  of  M-ary  FSK  and  QPSK  modulation  are  combined,  a  hybrid  digital 
modulation  scheme  results  where  signals  are  both  phase  modulated  and  frequency 
shifted.  This  new  modulation  technique  will  have  four  signals  at  each  of  the  operating 
frequencies  corresponding  to  the  four  different  possible  phase  modulations.  The  signal 
set  can  be  expressed  as 

4*)  =  ypf  cos[2jt/,/  +  /y]  (3.1) 

for  i  =  0,  1,  2,  3  and  0  <;  r  <  Ts ,  and 

4')  =  cos[27r/2/  +  (/  -  4)  f  ]  (3.2) 

for  i  =  4,  5,  6,  7  and  0  <,  t  £  T, . 

Integrating  the  square  of  the  signal  ^(r)  over  the  interval  0  <i<T,  yields 


=  |  ~cosW,/]rf/ 

-  J",  5in^4"^  ' 

L  w. 


(3.3) 


If  the  product  fT,  is  an  integer,  or  if/jr,  >  >  1  ,  then 


*T, 

*oW  =  E  (3-4) 

*4) 

which  corresponds  to  the  signal  energy.  Similar  integration  yields  identical  results  for 
s,(/)  ,  i  =  1,  2,  ...,  7. 
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Applying  the  Gram-Schmidt  orthonormalization  procedure  Ref.  5  (pp.  266-269]  to 
the  signal  set  described  by  Eqs.  3.1  and  3.2  generates  a  set  of  orthonormal  basis  func¬ 
tions  having  the  property 


i;  for  '=/! 

0;  for  i 


(3.5) 


so  that  the  signals  s,{t)  can  be  synthesized  using  a  linear  combination  of  basis  functions, 
i.e., 


;=i 


where  the  s,,  coefficients  are  defined  as 


N 

'  =  °> 1 . 7 


(3-6) 


•'A 


(3.7) 


The  first  basis  function  is  obtained  by  normalizing  s„(r) .  From  Eq.  3.4,  we  have 


«£,(/)  =  -^S-  =  /-f  cos[27r/,/] 


(3.8) 


so  that  s0(t)  =  v'£  0](/),  and  the  coefficient  s0j  =  v  £  .  From  Eq.  3.7,  s„  is  given  by 

in  =  f 

=  cos(2nft  + -j  cos{2nft)dt 

s 

=  J  ^  cost  y  ]  +  cos[27r(2/1)t  +  y  ]}rfr  =  0 


(3-9) 


Defining  the  auxiliary  function  y,(r)  as 


yfr)  =  sl(r)-s]]4>](r) 
= 


(3.10) 
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where  the  second  equality  follows  from  Eq.  3.9,  the  second  basis  function  can 
tained  by  normalizing  the  auxiliary  function.  That  is 


,  m  Vito  5i(') 

sit 


v£ 


yi 


~f~  cos(2*/,/  +  y ) 


Using  next  the  coefficients  s21  and  s22  to  obtain  the  auxiliary  function  y2(r)  as 

y2(r)  =  s2(t)  -  i210,(r)  -  s22<f)2(t) 


where 


521 


=  f 

Jo 

=  £  cos(2n r/,r  +  n)yJ^-  cos(2« fyt)di 

i'  e  r  ^ 

=  --  '{  cos(>]  +  cosO(2/;)/  +  n\}dt  = -%fE 

1  S  Jq 


and 


s22  =  f 

Jo 

=  cos(27r/,/  +  co%{2nf^t  +  ~  )dl 

0 


]  +  cos[2«(y,  )M--y-]}<*«0 


results  in  Eq.  3.12  becoming  equal  to 

y2(0  -  s2(r)  1- J E  <t>i{t) 

=  cos(2nf\l  +  n)  +  cos(2rt/tr) 


2£ 


[  -  cos(2jr/,r)  +  cos(27r/]t)]  =  0 


This  indicates  that  s2(/)  is  completely  described  as 


be  ob- 

(3-11) 


(3-12) 


(3.13) 


(3.14) 


(3.15) 
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s2(/)  =  -V£0,« 


(3.16) 


and  thus  no  new  basis  function  is  introduced  by  s2(r)  . 
The  coefficients  for  s3(r)  are  given  by 


*31 


m 


Jo 

=  +  cos(2nfxt)dt 

Jo  s 

=  ^Y~  I  {  cos[  -y-  ]  +  cos[2^(2f,)/  +  ~  ]}Jt  =  0 
s  Jn  " 


(3.17) 


and 


532  =  S^i)4>i{i)di 

Jo 

“  cos(2”y;' + ]\f¥,  cos(2"-/;' + f  k' 

VT  fT  _ 

=  —  {  cosO]  +  cos[>(2./;  )/]}<//  -  -  V£ 

1  s  *>0 

so  that  the  auxiliary  function  y3(r)  becomes 


>'3(0  =  s3(0  ~  S32<t>2(0 

=  ^J2^  cos(2nfi:  +  3f)+yJl£-  cos(2rt/,r  +  — ) 
=  cos{2nf\t +  dL  )  +  cos(2»yjf  +  y)]  =  0 


(3.18) 


(3.19) 


Thus,  as  with  s2(r),  s3(r)  introduces  no  new  basis  functions.  Continuing  this  procedure 
for  the  signal  s«(r)  ,  we  have 
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*41 


= 

— oo 

=  yj^T~  COS(2^)7f  C0S(2 nfxt)dt 

2-Je  rn 

=  — = —  cos(2nf2t)  cos{2nfxt)dl  =  0 

Jq 


(3.20) 


where  the  last  equality  is  due  to  the  assumption  that  the  frequencies/  and/  are  spread 
far  enough  apart  to  maintain  orthogonality  as  in  binary  FSK  modulation.  Similarly 


*42  —  f  *4(002(0^ 

J0 

-  yj^  cos(2rr/20>yy-  cos(2tt//  +  y )dt 
2  ~E  f 

=  —7 —  cos(27r/2/)  cos(2rr/t  +  >/r  =  0 


(3.21) 


so  that  the  auxiliary  function  y«(r)  becomes 


4 

)'4(0  =  *4(0  -  Yjsij<t>j(t)  =  s4{t) 


(3.22) 


v=i 

Normalizing  y4(r)  in  order  to  obtain  the  third  basis  function,  results  in 

Si(r) 


sar)  /  -> 

<t> 3(0  =  — =-  =  yl- y-  cos(27r/2r) 

v  ^  V  1 J 


(3.23) 


Next,  evaluating  sSj  for  j  =  1,  2,  3  yields 


*si 


=  f\(/)0,(rMr 
Jo 

=  Vf  cos(2rr/2r  +  y  )^T  cos{2nf]t)dt 
2  /£  r  7; 

=  -7—  ,cos(2tt/2/  +  y )  cos(27r/  t)dt  =  0 

/  j  Jn  Z 


(3.24) 


24 


and 


*52  =  !s5(t)4>2(t)dt 

Jo 

Vf  cos(2 nf2t  +  y  )yf£  cos(2 nfxt  +  y  )rfr  (3.25) 

=  — y-  j*  ‘cos{2nf2t  +  y )  cos(27r/,f  +  y  )<*  =  0 

where  again,  the  frequencies  yj  andf2  are  chosen  to  maintain  orthogonality  as  explained 
above.  Furthermore, 


=  f  !s5(i)4>i(t)dt 

Jo 

=  7fcos(2^  +  f)7XCos(2^ 

•'O  5 

-  J^l  cos[  ]  +  cos[27t(2/2)r  +  f  ]}Jt  =  0 


•553 


(3.26) 


so  that  with  all  three  of  the  coefficients  equal  to  zero,  s5(r)  must  be  orthogonal  to  all  of 
the  basis  functions  generated  thus  far.  The  fourth  basis  function  is  then  obtained  by 
normalizing  s5(/)  ,  that  is 


,  ,  s  *s(0 

04(O  =  -=-  = 

v  E 


~  cos(27r/2r  +  y ) 


(3.27) 


It  can  be  easily  shown  that  the  remaining  signals  s6(t)  and  s-{t)  can  be  expressed  as 

—  (3.28) 

s,M  «  -  -JE  4,,U) 

Thus  the  coefficients  s63  and  s-*  are  both  equal  to  -  J~E.  Since  the  set  {<£//)}  represents 
an  orthonormal  basis  in  a  vector  space,  the  signals  s,(t)  can  now  be  interpreted  geomet¬ 
rically  as  vectors  having  components  equal  to  the  s,y,  namely 

5,(0  -*•  St  =  (sa,  5(2 . 5i  V) ;  i  =  0,  1,  . . . ,  M  -  1  (3.29) 
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where  N  is  the  number  of  linearly  independent  orthormal  basis  vectors,  and  M  is  the 
number  of  signal  vectors.l  It  is  well  known  that 


W)2-  l*l2  =  £, 


(3.30) 


where  E,  is  the  energy  of  the  ith  signal.  For  the  case  of  2FSK/QPSK,  M  is  equal  to  8, 
N  is  equal  to  4,  and  the  vectors  representing  the  signals  are 


so  =  (  +  JE  ,_0,  0,  0) 

s,  =  (0,  +JE,  0,  0) 

52  =  (  -  v'£  ,  0,  0,  0) 

s3  —  (0.  —  JE ,  0,  0) 

2  ’  (3.31) 

s4  =  (0,  0,  -f-  v  E  ,  0) 

s5  =  (0,  0,  0,  +  S'E ) 

s6  =  (0,  0,  -  %/£ ,  0) 

s-  =  (0,  0,  0,  -v/£) 

Since  there  are  eight  different  possible  transmitted  signals,  each  signal  can  represent 
three  bits  of  information. 

The  average  energy  and  also  the  dimensionality  of  the  signal  set  can  be  minimized, 
by  subtracting  from  each  signal  a  constant  vector  a  such  that 

.v-i 

2]^}|5;-a|2  (3.32) 

/= o 

is  minimized,  where  P(s,}  is  the  a  priori  probability  than  the  signal  s,(r)  has  been  trans¬ 
mitted.  Such  vector  a  can  be  obtained  from 

,Vf-l 

«  =  £/{*}} ;  =  E(?}  (3.33) 

1=0 


1  The  general  practice  of  labeling  the  number  of  signals  as  M  is  followed  here,  but  the  letter 
M  in  "MI:SK,QPSK"  represents  the  number  of  frequency  shifts. 
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For  a  signal  set  of  the  type  described  above,  assuming  all  a  priori  probabilities  P(s,}  are 
equal,  Eq.  3.33  for  a  reduces  to 


a 


M-\ 


m 


0 


(3.34) 


so  that  there  is  no  minimum  energy  transformation  that  would  reduce  the  dimensionality 
or  the  average  energy  of  the  signal  set. 

Figure  10  on  page  28  shows  the  vector  representation  of  the  2FSK/QPSK  signal 
set  shown  as  eight  vectors  in  a  four  dimensional  vector  space.  A  signal  set  consisting 
of  X  orthogonal  signals  augmented  with  the  negative  of  each  signal  is  referred  to  as  a 
biorthogonal  signal  set.  As  shown  by  Figure  10,  the  2FSK/QPSK  signals  form  this  type 
of  set.  For  a  general  biorthogonal  signal  set,  the  number  of  dimensions  of  the  vector 
space  is  always  equal  to  half  the  number  of  signal  vectors.  That  is 

M  =  2;Y  (3.35) 

while  the  number  k  of  bits  represented  by  each  signal  is  given  by 

k  =  log  2\I  (3.36) 

Finally,  the  number  Kf  of  signal  frequencies  in  the  MFSK  QPSK  scheme  is  one-fourth 
the  number  of  signals  or 


M  —  4Kf  (3.37) 

The  best  assignment  of  bits  to  the  M  signals  of  the  biorthogonal  signal  set  is  to 
match  those  signals  separated  by  the  greatest  Euclidean  distance  in  the  vector  space, 
with  symbols  having  the  greatest  Hamming  distance.  With  this  rule  in  mind  the  bit  as¬ 
signments  become  obvious.  That  is,  signals  which  are  antipodal  are  assigned  symbols 
which  correspond  to  binary  complements.  Once  this  assignment  is  done,  all  other  vec¬ 
tors'  Euclidean  distances  are  equal,  so  that  any  other  assignment  patterns  can  be  freely- 
chosen.  A  typical  bit  assignment  is  shown  for  the  case  of  M  =  8  in  Figure  10. 

B.  OPTIMUM  RECEIVER 

With  the  basis  functions  having  been  specified,  available  analytical  results  on  the 
structure  and  performance  of  receivers  for  discriminating  biorthogonal  signals  received 
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in  additive  white  Gaussian  noise  make  it  possible  to  design  an  optimum  receiver  for  the 
discrimination  of  MFSK'QPSK  signals.  The  received  signal  r(t)  is  assumed  to  consist 
of  a  signal  component  s,(t)  plus  a  sample  function  of  a  white  Gaussian  noise  process  n(t). 
That  is, 


r(c)  =  Si(t)  +  n(t)  (3.38) 

As  shown  by  Ref.  5  [pp.  211-235]  the  optimum  receiver  takes  on  the  structure  of  a  cor¬ 
relation  receiver.  A  block  diagram  of  the  general  correlation  receiver  is  shown  in 
Figure  1 1  on  page  30.  The  signal  r(t)  is  processed  by  a  bank  of  correlators,  where  the 
number  of  correlators  in  the  receiver  is  equal  to  N,  the  number  of  dimensions  in  the 
signal  vector  space  model.  The  effect  of  the  additive  white  Gaussian  noise  is  to  add  to 
the  signal  vector  a  noise  vector  n  having  N  independent  components,  each  of  which  is 
a  zero  mean,  Gaussian  random  variable  with  variance  .  Thus  the  correlator  outputs 
tj  form  the  components  of  a  random  vector  r,  i.e., 

r  =  s,  +  n  *  (r„  r2,  . . . ,  rv)  (3.39) 

After  generating  the  vector  r  ,  the  components  r;  are  processed  through  a  weighting 
matrix  to  generate  the  dot  products  d,  which  are  given  by 

v 

=  r  '%  =  Yjrjs‘J  (3-40) 

j=  i 

The  d,  are  then  summed  with  bias  terms  c,  where 

ct  =  y  (Ao  ln[P{^}]  -  1 %  1 2) ;  i  =  0,  1 . M  -  1  (3.41) 

Finally,  all  summer  outputs  1,  are  routed  to  circuitry  which  decides  which  signal  was  sent 
according  to  the  rule: 


If  !,>!*;  V  *  =  0,  1,  ...,  M—  1  ;  k  #  / 

(3.42) 

Then  decide  Si{t)  was  the  transmitted  signal. 

For  the  signal  vectors  of  Eq.  3.31,  assuming  P{s,(/)}  for  all  /,  the  structure  of 
the  optimum  receiver  simplifies  from  that  of  the  general  correlation  receiver  to  the 
2FSK  QPSK  receiver  shown  in  Figure  12  on  page  32.  The  2FSK  QPSK  scheme  can  be 
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extended  beyond  the  case  considered  by  increasing  the  number  of  signal  frequencies. 
For  example,  in  order  to  increase  the  number  of  bits  per  signal  from  three  to  four,  the 
number  of  signals  M  must  be  doubled  to  16.  This  is  accomplished  by  doubling  the 
number  of  frequencies  Kf  from  two  to  four,  and  thus  doubling  the  number  of  basis 
functions  from  four  to  eight.  This  change  will  be  reflected  in  the  receiver  structure  by 
a  doubling  of  the  number  of  correlators  to  eight.  The  new  signal  set  is  also 
biorthogonal,  and  further  extensions  to  MFSK/QPSK  result  in  N  basis  functions,  Kf 
frequencies,  M  signals,  and  k  bits  per  signal.  The  formula  for  the  probability  of  a  correct 
symbol  decision  for  the  biorthogonal  set  is  given  by  Ref.  5  [pp.  261-263]  as 


r  oo 

PsiQ = J7n(«  -  -Jew 1 1  -  2j7.ow 


(3.43) 


where  fn(x)  is  a  Gaussian  distribution  function  expressed  as 


/nM  = 


and  where 


(3.44) 


(3.45) 


a-y'E 

Using  the  transformation  y  = - 5 - and  rearranging  the  terms  in  Eq.  3.43  yields 


PsiQ  = 

*oc 

1  —  2  [°°  Mdp 

% 

Jay+  yJ~E 

-i-v-i 


g(y)dy 


Using  another  variable  transformation  the  inner  integral  becomes 

J*  OO  Poo  /  / 

MP)dp  =  g{x)dx  =  (2)’  +  ^ 

-v+./F  J ...  VE  V 


'<ry+  y/T 

Substituting  this  result  into  Eq.  3.46  gives 


E_ 

a 


(3.46) 


(3.47) 
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Figure  12.  Diagram  of  Correlation  Receiver  for  2FSK/QPSK 
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(3.48) 


Thus  we  have  the  probability  of  a  correct  symbol  decision  associated  with  the  optimum 
receiver  expressed  in  terms  of  the  Q  function. 

C.  OPTIMUM  RECEIVER  BER 

In  order  to  compute  the  Bit  Error  Rate  (BER)  of  the  optimum  receiver  shown  in 
Figure  12,  it  is  necessary  to  determine  under  what  conditions  a  chosen  bit  will  be  cor¬ 
rectly  indentified  by  the  receiver  for  a  given  transmitted  signal.  The  three  bits  assigned 
to  each  2FSK  QPSK  signal  are  identified  as  IB,  2B  and  3B  corresponding  to  the  most 
significant  bit,  the  second-most  significant  bit,  and  the  least  significant  bit  respectively. 
We  can  assume  without  any  loss  in  generality  that  sQ(t)  is  transmitted.  The  probability 
that  the  IB  is  correctly  identified  at  the  receiver,  conditioned  on  the  assumption  that 
s0(f)  was  transmitted  can  be  written  as 


Pr{\B  Correa |70}  =  /V{r,  =  max  r|70}  4-  Pr{  -r2  -  max  r,|70} 

j  J  _  j  J 

+  Pr{ r3  =  max  r  I  s0}  +  Pr{  -r4  =  max  r  l?0} 

J  j 


(3.49) 


Now 

Pr{ r,  =  max  ry|70} 

j  _  (3.50) 

■  Pr\r\  >  -rj.  r,  >  r2,  r,  >  -r2,  r,  >  r3,  r,  >  -r3,  r,  >  r4,  r,  >  -r4|s0} 

If  we  further  condition  on  r„  that  is  r,  =  R,  ,  since  r,  >  -r,  only  for  R{  >  0  ,  then  this 
conditioning  will  be 

=  max  r,-|s0,  r,  =  Rj  >  0} 

j  _  (3.51) 

=  Rr{R,  >  r2,  R,  >  -r2,  R,  >  r3,  R,  >  -r3,  R,  >  r4,  R,  >  -r4|s0,  r,  =  R,  >  0} 

Observe  that 

Rj  >  Tj,  R,  >  -Tj  =>  -R,  <  Xj  <  R,  (3.52) 


so  that 
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Pr{  -/?,  <  Xj  <  Rl  |  s0,  r,  =  R}  >  0}  = 


=  2 


V 

_1_ 

2 


/  2na2 


exp 


z5. 

la2 


\dR, 


f*OC 


C 


V 


■p="  exp(  — ~  )dx 
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Pr{  r,  =max  r,|70} 
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3^rexp(^  » 


(3.59) 


=  l-2g(4L);  for  j  =  2,  3,  4 
Using  this  result  with  Eqs.  3.53  and  3.58  yields 

/V{r,  =  max  r,|s0,  r,  =  R]  >  0}  =  ^  1  -  2g(  )]  (3.60) 

so  that  integrating  over  all  Rl  >  0  yields 


(3.61) 


Note  that  this  expression  approaches  unity  as  —  -» oo  as  expected.  Consider  now 
Pr{  — r2  =  max  r )  I  s0} 

'  (3.62) 

=  Pr{  — r2  >  r, ,  -r2>-r„  -r2  >  r2,  -r2  >  r3,  -r2  >  -r3,  -r2  >  r4,  -r2>-r4|s0} 


If  we  further  condition  on  r2,  that  is  r2  —  R2  ,  since  —  r2  >  r2  can  only  be  satisfied  for 
R2  <  0  ,  then  this  conditioning  will  be 


Pr{  -r2  =  max  r,  1 s0.  R-i<  0} 

j  ^  (3-63) 

=  Pr{ r,  <  -R2,  r,  >  R2,  r3  <  -R2,  r3  >  R2,  r4  <  -R2,  r4  >  R2  |s0,  r2  =  R2  <  0} 


Let  R'2  =  —  R2 ,  so  that  in  terms  of  R'2  , 
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rj<-R2,rj>R2  =>  rj<R'2,  rj>-R’2  =*  -R'2  <  ry  <  R'2  ;  j=  1, 
and  the  requirement  R2<  0  =>  R'3>0  ,  so 

Pr{  -r2  =  max  ry|  s0,  r2  =  -R'2,  R'2  >  0} 

4 

=  -R'2  <  r )  <  /?'2 I*,,  r2  =  -/?'2.  R'2  >  0} 

M  i 
M2 

Foij=  1  , 


/V{  -tf'2  <  r,  <  /?'2} 

ru'i 
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since  /?'2  >  0  .  For,/  =3,4 


so  that 


-/?'2  <  ry  <  7?’2}  = 


-a; 


UV2™2 


exp 
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2  /^y 


°  .  /  2 
— p=-exp(  — r—  Idx 
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3,  4(3.64) 


(3.65) 


(3.66) 

\dx 


(3.67) 
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Pr{  -r2  =  max  r;  |T0,  r2  =  -R'2,  R'2  >  0} 


=  ;,-22(-^)] 


l-Q 


R'2  +  JE 

a 


-Q 


R'2-JE 


(3.68) 


Since  r2  takes  on  values  R2,  we  can  define  a  new  random  variable  r'2  =  -r2  where  r'2  takes 
on  values  R'2,  i.e.,  if  r2  =  R2,  then  r'2  =  —  R2=  R'2 .  So  the  condition 

r2  =  -R'2,  R’2>  0  =>  -r'2  =  -R’2,  R'2>  0  =>  r'2  =  R'2,  R'2  >  0  (3.69) 

and  given  the  symmetric  p.d.f.  for  r2  conditioned  on  s0 ,  we  see  that  the  conditional  p.d.f. 
of  r'2  is  no  different  than  the  conditional  p.d.f.  of  r2 ,  so 


Pr{  — r2  =  max  t)-l70}  =  Pr{ r’2  =  max  ryls0) 


['-22(4)] 


1-2 


R'2  +  JE  \  J  R’2-JE 


-Q 


,  (3.70) 

\fr(R'2)dR'7 


Through  a  transformation  of  variables  this  expression  is  reduced  to 


/*r{  — r2  =  max  r,|70} 

J 

foo  2 

=  d -2  20): 


1  /VW-7" 


J  v 


==-  exp(  )dy 
1 2n 


Note  that  this  last  expression  approaches  zero  as  —  -+  oo  .  Consider  now 


/V{r3  =  max  rj\s0} 

=  Pr{ r3  >  r„  r3  >  -rt,  r3  >  r2,  r3  >  -r2,  r3  >  -r3,  r3  >  r4,  r3  >  -r4 1 (3-72) 
-  <  r3>  ri  >  ~r3-  r2  <  r3»  r2  >  — r3,  2r3  >0,  r4  <  r3,  r4  >  -r3 I 

If  we  condition  on  r3  that  is  r3=  R%  ,  we  must  have  f?3>0,  or  we  cannot  satisfy  the 
2r3  >  0  condition.  Thus 


Pr{ r3  =  max  r,  I  s0,  r3  =  R3  >  0} 
j 

=  Pr{  —R3  <  f]  <  R3,  —R3  <t2<  R3,  —R3  <  r4  <  R3 1  s0.  r3  —  ^3  > 


(3.73) 


37 


From  the  mathematical  form  of  this  equation,  it  is  obvious  that  this  probability,  when 
averaged  over  all  /?3  >  0  produces  the  same  probability  as  that  given  by  Eq.  3.71,  which 
also  is  identical  to  Pr{  — r4  =  max  ry|s0}  .  Thus,  we  obtain  from  Eqs.  3.49,  3.61  and  3.71 


Pr{  1 B  Correct  |  s0} 

r°°  / —  3 

=  I  |^1  -2Q(y  +  — g(y)dy 


i 


& 

<7 

OO 


(3.74) 


+  3  [1-2  <20011 


i-e(y+ 
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g{y)dy 


If  the  2B  is  selected  then 

Pr{2B  Correct  I  s0}  =  Pr{ rt  =  max  rylsg)  4-  /V(r2  =  max  iy|s0} 

+  Pr{ r3  =  max  r,  I  s0}  +  Pr{ r4  =  max  r< 1 

j  j 

The  first  term  on  the  right  hand  side  of  Eq.  3.75  is  given  by  Eq.  3.61,  while  the  last  three 
terms  are  all  of  the  same  form  as  Pr{ r3  =  max  r,|s0}  from  Eq.  3.72  and  will  all  yield  the 
same  probability  as  given  by  Eq.  3.71.  Therefore,  the  probability  that  the  2B  is  correct 
is  the  identical  to  the  probability  that  the  IB  is  correct.  It  can  be  shown  in  similar 
fashion  that  the  probability  of  correctly  identifying  the  3B  is  also  equal  to  that  for  the 
IB.  Thus  Eq.  3.74  is  also  the  expression  for  Pb{Conect\ s0)  the  average  probability  of 
decoding  a  bit  correctly  given  that  s0(t)  was  transmitted. 

Now  consider  the  case  of  a  given  signal  sk{t)  being  transmitted,  where  sk(t)  can  be 
any  one  of  the  signals  s,(r),  /  =  0,  1,  . . . ,  7  .  For  any  chosen  bit,  the  probability  that 
this  bit  is  correct  conditioned  on  the  hypothesis  that  sk(t)  was  transmitted  can  be  ex¬ 
pressed  in  a  form 


4 

Pr{ Chosen  Bit  Correct 1 7k}  =  Pr{ r,  =  max  Tj\7k}  +  ^Pr{rm  =  max  tj 1 7k)  (3.76) 

J  m=l  J 

m*l 

where  r,  and  the  r„  terms  are  the  signed  correlator  outputs  ±  r; .  The  output  r,  corre¬ 
sponds  to  sk  in  that  if  r,  is  the  maximum  of  all  the  ±  r;  ,  then  the  receiver  decides  sk{t ) 
was  transmitted.  The  outputs  rm  correspond  to  three  other  s<  which  if  selected  by  the 
receiver,  would  result  in  the  chosen  bit  still  being  correct.  The  first  term  on  the  right 


38 


hand  side  will  again  yield  the  same  probability  as  that  given  by  Eq.  3.61.  The  terms  in 
the  summation  are  all  of  the  same  mathematical  form  as  Eq.  3.72,  and  so  are  equal  to 
the  probability  given  by  Eq.  3.71.  Thus  Eq.  3.74  is  the  expression  for  the  probability 
that  any  chosen  bit  is  correct  for  any  given  signal  transmitted,  and  therefore  yields  the 
overall  average  probability  of  a  correct  bit,  Pt{C)  . 

Some  insight  into  the  total  probability  of  a  correct  bit  decision  P4{C),  may  be  gained 
by  comparing  Eq.  3.74  with  the  probability  of  a  correct  symbol  decision  Pt{C}  given  by 
Eq.  3.48.  Since  N,  the  number  of  basis  functions  is  equal  to  4,  the  first  term  in  the  ex¬ 
pression  for  Pb{C }  is  equal  to  P,{C}.  This  should  come  as  no  surprise,  since  the  deriva¬ 
tion  for  the  probability  of  the  IB  being  correct  began  with  Eq.  3.43  where  the  first  term 
is  in  fact  the  probababilty  that  this  bit  is  correct  because  the  receiver  made  the  correct 
symbol  decision. 

The  second  term  in  Eq.  3.74  represents  the  probability  that  the  bit  is  correctly 
identified  even  though  the  receiver  has  made  a  symbol  error.  Recall  from  Figure  10  that 
only  those  signals  which  have  vector  representations  orthogonal  to  each  other  have  bits 
in  common,  and  for  2FSK,QPSK  there  are  three  such  signals  meeting  this  requirement 
for  any  one  fixed  signal.  The  probabilities  of  erroneously  choosing  one  of  these  signals 
over  say  s0(/)  ,  are  all  equal  as  discussed  above.  Furthermore,  the  conditional  probability 
of  (incorrectly)  choosing  any  of  the  signals  orthogonal  to  s0(/),  can  be  shown  to  have  the 
mathematical  form  of  Eq.  3.72.  Thus  Eq.  3.71  specifies  the  probability  of  erroneously 
choosing  any  one  of  the  signals  orthogonal  to  j0(r) . 

The  above  results  for  2FSK  QPSK  can  be  extended  to  \1FSK  QPSK,  that  is,  a 
multi-frequency  four- phase  signaling  scheme.  The  structure  of  the  optimum 
MFSK  QPSK  receiver  is  the  same  as  that  of  the  general  correlation  receiver  shown  in 
Figure  11  on  page  30  where  the  bank  of  N  correlators  corresponds  to  a  pair  of 
correlators,  one  each  for  the  In-phase  and  Quadrature  channels,  for  each  of  the  K,  fre¬ 
quencies.  The  expression  for  the  conditional  probability  of  a  correct  signal  decision  has 
already  been  given  in  Eq.  3.48.  The  expression  for  the  conditional  probability  of  decid¬ 
ing  in  favor  of  a  signal  that  is  orthogonal  to  the  transmitted  signal,  yet  does  not  result 
in  a  bit  error  takes  on  a  form  similar  to  that  derived  for  the  2FSK/QPSK  case.  To 
demonstrate  this,  assume  that  s0(r)  has  been  transmitted  and  examine  the  probability 
that  one  of  the  signals  orthogonal  to  s0(r)  is  erroneously  choosen  over  s0(r ) .  This  corre¬ 
sponds  to  the  probability  that  for  some  r* ,  k  =  2,  3,  ... ,  N 
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Pr{rk  =  max  ry|s0} 

=  Pr{xk  >  rlt  rk>  -r„  rk  >  r2,  rk  >  -r2, 


(3-77) 

rk  >  -rk,  rk>  r,v,  r*  >  -r,v|s0} 


Further  conditioning  on  r*  yields 


Pr{rk  =  max  ry  1 70,  rk  =  Rk,  Rk  >  0} 
v 

= r>{ <  Rk  i  *  ** = **  >  o}  (378) 

7=1 

Mk 

The  First  term  in  the  summation  (/  =  1)  is  the  same  as  that  given  by  Eq.  3.66  for 
2FSK  QPSK.  All  the  other  terms  (j  =  2,  3,  . . . ,  A';  j  ¥=  k)  are  given  by  Eq.  3.67,  so 
that 


Pr{rk  =  max  r y  I  s0,  r*  =  Rk,  Rk  >  0} 


-['-22(4L) 
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(3.79) 


Now,  integrating  over  all  r*  >  0  yields 


Pr{rk  =  max  r,|70} 
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=  Cl  -2<20)] 


-V-2 
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Since  there  are  M  —  2  signals  orthogonal  to  s0(r)  ,  and  for  any  chosen  bit,  half  of  these 
signals  have  this  chosen  bit  in  common  with  the  one  in  .%(/),  yields  ^ ^  ^  =  N  —  1 
ways  of  making  a  symbol  error  while  correctly  indentifying  the  bit  in  question.  Thus  the 
expression  for  the  probability  of  a  correct  bit  decision  conditioned  on  the  signal  s„(/) 
being  transmitted,  is  given  by 
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Using  steps  similar  to  those  used  in  the  analysis  of  the  2FSK/QPSK  receiver  as  well  as 
those  presented  above,  it  can  be  shown  that  Hq.  3.81  is  also  the  expression  for  the  overall 
average  probability  of  a  bit  being  correctly  identified  by  the  MFSK/QPSK  optimum  re¬ 
ceiver  Pt{C}  . 

Since  the  integrals  involved  in  calculating  the  bit  error  probability  for  the 
MFSK  QPSK  optimum  receiver  cannot  easily  be  expressed  in  closed  form,  the  BER  can 
only  be  evaluated  using  numerical  techniques.  Appendix  B  presents  a  computer  pro¬ 
gram  in  FORTRAN  language  used  to  compute  the  actual  bit  error  probability  for  the 
optimum  receiver  by  employing  numerical  integration.  An  upper  bound  for  the  receiv¬ 
er's  bit  error  probability  is  given  by  Ref.  6  [pp.  255-257]  as  a  union  bound,  namely 


Pb{e}<Y 


(3.82) 


This  bound  is  a  good  approximation  to  the  actual  BER  for  k  >  3  at  moderate  to  laree 

A 

A  o ’ 

D.  DBD  RECEIVER  FOR  MFSK/QPSK 

Another  receiver  structure  for  MFSK  QPSK  signals  is  implementable  by  using  DBD 
techniques.  Ref.  2  [pp.  19-25]  A  diagram  of  the  structure  of  this  receiver  is  shown  in 
Figure  13  on  page  42  for  the  2FSK/QPSK  case.  The  DBD  receiver  uses  the  same 
correlator  subsystem  as  the  optimum  receiver  Figure  12  on  page  32,  however  the  two 
receivers  differ  in  the  way  the  correlator  outputs  are  processed.  Whereas  the  optimum 
receiver  processed  these  correlator  outputs  through  a  weighting  network  and  then  se¬ 
lected  its  largest  output,  the  DBD  receiver  routes  the  correlator  outputs  through  a  two 
stage  summing  network  and  then  compares  the  summer  outputs  to  a  zero  threshold  in 
order  to  arrive  at  binary  decisions. 
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Figure  13.  Diagram  of  DBD  Receiver  for  2FSK/QPSK 


E.  MFSK/QPSK  DBD  RECEIVER  BER 

The  analysis  for  the  derivation  of  the  BER  performance  for  the  2FSK/QPSK  DBD 
receiver  is  similar  to  that  carried  out  for  the  performance  analysis  of  the  optimum 
(symbol  discrimination)  receiver.  We  begin  by  assuming  ^(r)  is  transmitted.  Since  the 
correlator  outputs  are  identical  for  both  receivers  in  question,  the  analysis  on  the  per¬ 
formance  of  the  DBD  receiver  can  proceed  with  the  correlator  outputs,  which  are  again 
characterized  as  independent  Gaussian  random  variables  with  density  functions  given 
by  Eqs.  3.56  and  3.57.  Note  from  Figure  13  that  there  is  a  path  from  each  correlator 
output  to  the  input  of  all  second  stage  summers.  The  outputs  of  these  summers,  denoted 
b„ ,  n  =  1,  2,  3  ,  are  given  by 


b,  =  r!  -  r2  +  r3  —  r4 

b2  -  rj  +  r2  +  r3  +  r4  (3.83) 

b3  =  ~ri  -  r2  +  h  +  r4 

Each  b„  will  also  be  a  Gaussian  random  variable  having  a  mean  equal  to  the  sum  of  the 
means  of  the  r„  and  due  to  the  independence  of  the  r,,  the  variance  of  b„  will  consist  of 

jY 

the  sum  of  the  variances  of  the  r;  ,  each  of  which  equals  — Thus  we  have 

b,  ~  A V£ ,  2A0) 

b2  ~  ,V(v£  ._2A0)  (3.84) 

b3  ~  A'(  —  v  £  ,  2V0) 

where  N{a,  b)  is  used  to  denote  a  Guassian  probability  density  of  mean  a  and  variance 
b.  With  the  statistics  of  the  b„  specified  comparing  these  random  variables  to  a  zero 
threshold  yields  the  following  conditional  probabilities  of  decoding  the  bits  correctly 
namely 


Pr{15  Correct]  s0)  ~  Pr(b5  >  0}  = 


*oo 

(x-JE)2  ' 

J„  JM2\)  xp 

2(2Ao) 

dx 


(3.85) 
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(3.86) 


Pr{2B  Correal  s0]  =  Pr{ b2  >  0}  = 


/*oo 

1  1 

{x-jEf  ' 

I  . .  —  CAD 

J0  v'2*(2Ao) 

2(2.V0) 

dx 


Pr{1B  Correct  |  s0)  =  /V{b3  <  0}  = 


=  1-0 


(x  +  JE)2  ' 

J  „'2*(2,V0)  p 

2(2  N0) 

dx 


(3.87) 


2A0 


Since  each  of  the  individual  bits  have  identical  probability  of  correct  decoding,  the  cor¬ 
responding  average  probability  conditioned  on  the  assumption  that  s0(t)  was  transmitted 
is 

Pb{Correa\%)  =  1  -  (3-88) 

and  consequently  the  conditional  bit  error  probability  is 

Pb{Error\70)  =  1  -  Pb{Correa |?0)  =  (3-89) 

Consider  next  the  case  in  which  any  one  of  the  s,(t) ,  i  =  0,  1,  ....  7  is  transmitted,  for 
any  of  the  receiver  correlators.  If  the  transmitted  signal  is  matched  to  the  local  signal 
correlator,  its  output  r,  will  be 

Tj  =  JE  +  n,  (3.90) 

and  for  the  correlator  having  antipodal  local  signal  matching,  the  output  will  be 

r,  =  -  J~E  +  tij  (3.91) 

All  other  correlators  will  have  an  output  given  by 

Xj  =  n;  (3.92) 
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Examining  now  any  one  of  the  second  stage  summer  outputs  b„ ,  if  the  bit  corresponding 
to  b„  for  the  signal  sent  is  a  logical  zero,  then  b„  is  statistically  characterized  by 


bn~AV£,  2.Y0)  (3.93) 

In  this  case  the  bit  will  be  correct  if  b„  is  greater  than  the  zero  threshold,  and  the  prob¬ 
ability  of  the  bit  corresponding  to  b„  being  correctly  decoded  is  given  by 


(3.94) 


However,  in  the  case  where  the  bit  is  a  logical  one,  then 

bn~.Y(-,/£.  2A'0)  (3.95) 

Now  the  bit  will  be  correct  if  b„  is  less  than  the  zero  threshold,  and  the  probability  of 
correct  decoding  is  given  by 


/V{b„<0} 


V2*(2.Y0) 


=  1  - 


(*  +  vE)2 

2(2-Y0) 


(3.96) 


For  both  cases  the  probability  of  a  bit  error  is  the  same,  and  it  is  identical  for  all  b„  re¬ 
gardless  of  which  signal  s,(r)  is  transmitted.  Thus,  we  can  write  the  total  bit  error  prob¬ 
ability  for  the  2FSK.  QPSK  DBD  receiver  as 

(M7) 

Extending  the  2FSK  QPSK  results  to  the  general  case  of  MFSK/QPSK  in  which 
there  are  M  signals,  N  receiver  correlators,  k  bits,  and  Kf  frequencies  where  the 
biorthogonal  signal  relationships  of  Eqs.  3.35  -  3.37  again  hold,  for  each  of  the  random 
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variables  b„ ,  n=  1,  2,...,  k  ,  the  corresponding  variance  will  be  the  sum  of  the  variances 
of  the  n j,  so 


V  Ao  A'A’n 

Var{bn)  =  =  ~T~  “  (3l98) 

>1 


and  again  the  mean  of  the  b„  will  be  ±  JE  .  Following  the  same  steps  as  those  carried 
out  in  obtaining  Eqs.  3.90  -  3.97,  the  total  bit  error  probability  for  MFSK/QPSK  DBD 
receiver  is  given  by 

(3"> 


F.  BANDWIDTH  OF  MFSK/QPSK 

The  bandwidth  of  MFSK  QPSK  signals  can  be  derived  from  knowledge  of  QPSK 
and  MFSK  bandwidth  requirements.  Each  carrier  is  switched  at  the  symbol  rate  Rs 
where  Rs  =  ~r  .  This  gives  each  carrier  a  two-sided  bandwidth  of  2 Rs  Hz.  To  maintain 

*  s 

orthogonality  between  signals  at  different  frequencies  in  MFSK  modulation,  the  fre¬ 
quencies  are  separated  by  Rs  Hz.  This  same  relationship  must  also  be  maintained  in 
MFSK  QPSK.  The  spectrum  of  the  2FSK  QPSK  signal  is  plotted  in  Figure  1-4  on  page 
47.  From  the  graph,  the  total  null-to-null  bandwidth  of  the  2FSK  QPSK  signal  can  be 
seen  to  be 


IF  =3  Rs  (3.100) 

Proceeding  to  the  family  of  MFSK/QPSK  signals,  Figure  15  on  page  48  shows  the 
spectrum  for  such  signals.  By  examining  this  figure,  the  total  null-to-null  bandwidth  can 
be  shown  to  be 


W'-^+QK, 

=  (2k~2  +  l)/?f 


(3.101) 


This  general  result  for  MFSK  QPSK  correctly  yields  the  result  for  the  special  case  of 
2FSK  QPSK  in  which  k  —  3,  namely 

tF  =  (23_2  +  !)/?,=  1RS  (3.102) 
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Figure  14.  Magnitude  Spectrum  for  2FSK/QPSK:  Dashed  lines  show  individual 

spectral  components. 


4S 


which  agrees  with  Eq.  3.100. 

The  bandwidth  efficiency  of  a  digital  modulation  scheme  is  given  by  Ref.  6  [pp. 
393-395]  as 


r i  = 


(3.103) 


and  is  one  way  to  evaluate  how  efficiently  different  modulation  schemes  use  the  available 
channel  bandwidth.  For  M-ary  PSK,  bandwidth  efficiency  is  given  by  Ref.  3  [pp. 
211-213] 


MPSK  2  R 

kRs 

=  2  Rs 
2 


I 


(3.104) 


and  for  M-ary  FSK  the  bandwidth  efficiency  is 

Vfsk  (M+l)Rs 
kRs 

~  (-V/+ 
k 

2k+  1 


(3.105) 


For  MFSK  QPSK,  using  Eq.  3.101  we  obtain  a  bandwidth  efficiency  given  by 


>1  MFSK! QPSK  '' 


(2k~2+Ws 

kRs 

(2k~2  +  1  )RS 


(3.106) 


2k~2  +  1 


A  table  of  bandwidth  efficiencies  for  M-ary  PSK,  M-ary  FSK,  and  MFSK, QPSK  is 
given  in  Chapter  5. 
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IV.  RAYLEIGH  FADING 


A.  THE  RAYLEIGH  FADING  MODEL 

Often,  when  a  signal  is  transmitted  through  a  channel  in  which  multipath  propa¬ 
gation  and/or  heavy  precipitation  conditions  exist,  fading  occurs  with  a  consequent 
change  in  the  signal  amplitude.  One  way  to  model  this  fading  is  to  redefme  the  signal 
amplitude  A  which  up  until  now  has  been  treated  as  a  constant,  as  a  Rayleigh  distributed 
random  variable  A  with  p.d.f. 

AM— !t«p( 

Ao  \  2A0  / 

where  u(-)  is  the  unit  step  function,  and  from  Ref.  7  [pp.  110-112],  the  value  of  Ac  is  re¬ 
lated  to  the  second  moment  of  A  as 

£{A2}  =  2Aq  (4.2) 

In  order  to  evaluate  the  effect  of  Rayleigh  fading  on  the  BER  for  a  given  modulation 
scheme  and  corresponding  receiver,  the  expression  for  BER  under  nonfading  channel 
conditions  must  be  averaged  over  the  distribution  of  the  signal  amplitude.  This  requires 
expressing  the  bit  error  probability  as  a  function  of  signal  amplitude,  weighting  the 
function  with  the  amplitude  p.d.f.  and  integrating  this  product  over  all  values  of  A. 

B.  EFFECT  ON  BER  OF  M-ARY  PSK  USING  DBD  RECEIVER 

For  the  S-PSK  DBD  receiver,  the  (non-fading)  bit  error  probabilities  are  given  by 
Eqs.  2.1  -  2.2.  These  probabilities  must  now  be  expressed  as  functions  of  the  random 
amplitude  A.  From  Eq.  1.1  it  is  apparent  that  the  signal  amplitude  is  given  by 

^  (43) 

so  that 

A2Ts 

E  =  -yi  (4.4) 

Since  the  signal  amplitude  is  now  a  random  variable,  the  signal  energy  also  becomes  a 
random  variable  E.  which  is  expressed  as 
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(4.5) 


Substituting  this  result  into  Eqs.  2.1  through  2.3  yields 
Pr{ IB  incorrect}  = 


i  cos  a) + 1  e(Avir sin  a 


Pr{ 2B  incorrect}  = 


q(^aJ-^[  cos  a  +  sin  a]^  +  cos  a  -  sin  a]^ 


and  as  before, 


2(21  A y -jy-  [  cos  a  +  sin  a]  j£>l  Ay  [  cos  a  -  sin 


Pr{3B  incorrect}  =  Pr{lB  incorrect} 


The  bit  error  probability  under  fading  conditions  for  the  most  significant  bit,  namely  IB, 
is  obtained  by  weighting  the  probability  given  in  Eq.  4.6  by  the  p.d.f.  of  the  random 
amplitude  A  and  integrating  over  all  A.  That  is 

Poo 

/^.{IB  incorrect}  =  I  P>{1B  incorrect ]fx{A)dA 

— QO 

00  i  (  r t  \  i  /  /  T  \ 

=  cos  a  J  +Y  Ay-y^  sin  a  )  fx{A)dA 

oot- 

This  integral  can  be  solved  using  the  identity  given  by  Ref.  8  [pp.  12-13]  as 


f  °°  Q{<*-Ayx{A)dA  =  y  !-• 

*  —  oo 


(4.10) 


Using  Eq.  4.2  and  the  linearity  of  the  expectation  operator  we  can  write  an  expression 
for  the  average  signal  energy  under  fading  conditions,  namely  £{E}  =  E  as 
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(4.11) 


=  -^-£{A!}=i(24) 
-  tsaI 


Using  Eq.  4.10  in  order  to  evaluate  Eq.  4.9  yields  the  bit  error  probability  under  Rayleigh 
fading,  namely 


Pxay{lB  incorrect}  = 

.  .r  .  i, 


(4.12) 


Next,  the  BER  under  fading  conditions  for  the  second-most  significant  bit,  namely  2B, 
is  obtained  by  integrating  the  product  of  the  BER  given  by  Eq.  4.7  with  the  Rayleigh 
p.d.fi.  The  contribution  from  the  first  two  terms  of  Eq.  4.7  can  be  derived  using  the 
same  steps  as  those  carried  out  above  for  the  evaluation  of  PRcj,{1B  incorrect}.  Since  the 
third  term  in  Eq.  4.7  involves  the  product  of  two  Q  functions,  a  new  identity  is  required. 
That  is 


r  oo 

QWA)Q{($AYx{A)dA  = 

^  —  OO 


1 

4 


/?  tan' 


2  n 


(4.13) 


the  proof  of  which  is  presented  in  Appendix  A.  Using  this  identity  for  the  contribution 
of  the  third  term  of  Eq.  4.7  and  combining  this  with  the  contributions  from  the  first  two 
terms  of  Eq.  4.7  yields 


PRay{2%  incorrect}  = 

1  1  1 


where  u{  and  u2  are  given  by 


u,  =  cos  a  +  sin  a 
u2  =  cos  a  —  sin  a 


(4.15) 


The  overall  average  bit  error  rate  under  fading  conditions  can  be  obtained  from 


pb,  i *ay{*}  =  y  ^(IB  incorrect)  +  y  PRa>,{2B  incorrect}  (4.16) 


For  the  16-PSK  DBD  receiver  operating  under  Rayleigh  fading  conditions,  much 
of  the  performance  analysis  can  be  carried  out  using  the  same  procedure  as  the  one  used 
for  the  evaluation  of  the  performance  of  the  8-PSK  DBD  receiver  operating  under  simi¬ 
lar  Rayleigh  fading  conditions.  The  bit  error  probabilities  under  non-fading  conditions 
are  given  by  Eqs.  2.4  -  2.6.  Rewriting  these  equations  in  terms  of  the  random  variable 
A,  for  the  most  significant  bit  we  obtain 


where  the  variables  v,  are  defined  as 


(4.18) 


v,  =  cos  a 
v2  =  cos  P 
v3  =  sin  a 
v4  =  sin  /? 

and  for  the  second-most  significant  bit 


/V{2B  incorrect}  = 


(4.19) 


where  the  variables  w,  are  given  by 


w,  =  cos  a  —  sin  a 
vv2  =  cos  a  +  sin  a 
w3  =  cos  p  —  sin  /? 
w4  =  cos  /?  +  sin  p 


(4.20) 


For  the  third-most  significant  bit 


/V{3B  incorrect}  = 

T 

A 


7' S' 

i=l 


A  T, 


2.V 


sfhr,  ot,  p)  Wr,  a,  p) 


(4.21) 


where  the  functions  /{T,  a,  p)  and  g,(r,  «,  p)  are  defined  in  Table  1  of  Chapter  2.  For 
the  least  significant  bit,  from  Eq.  2.3  it  is  obvious  that 

/V{4B  incorrect}  =  /V{1B  incorrect}  (4.22) 
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For  the  most  significant  bit,  the  BER  under  fading  conditions  is  obtained  using  Eq.  4.17 
and  the  identity  given  by  Eq.  4.9,  to  yield 


/^{IB  incorrect} 


-I-+Z 

i=i 


■ 


An 


=  F^^B  incorrect} 


-l 


+  1 


(4.23) 


where  the  variables  v,  are  defined  in  Eq.  4.18,  and  the  last  equality  in  Eq.  4.23  follows 
from  Eq.  4.22.  The  BER  under  Rayleigh  fading  conditions  for  the  second-most  signif¬ 
icant  bit  is  obtained  using  the  identities  of  Eqs.  4.9  and  4.13  to  yield 


PRay{ 2B  incorrect} 


4 

-ME 

i=  1 


E  2 


2.Y 


r~  W:  1  +  1 


iv,  tan 


-l 


+ 


Wi  +  2 


-i 


iv. 


w2  tan 


-l 


i  2  ,  J  E 

27T  /  w,  +  2[  — 


+ 


1  2  ,  ->/  E  ^ 
^  +  2\~, 


'0 


w. 


2 n  I  vv’2  -f  l{  ^ 


-l 


(4.24) 


where  the  w,  terms  are  defined  in  Eq.  4.20.  For  the  third-most  significant  bit,  the  inte¬ 
grals  to  be  evaluated  are  all  of  the  basic  form 
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M 

•'-OO  JC 


8expf  -4^/(r.  a,  P))gi(r,  a.  P)dTfx{A)dA 


-oo  •'O 

7T 


2-V0 


-  K/r,  «, 

**0  e/^o 

u 

=J>(r' 


x 2  exp(  —  x2)dxdT 


(4.25) 


where  c,  is  defined  as 


«/s  Pi  + 


2Ac 


(4.26) 


and 


a,  P) 


(4.27) 


The  variable  x  is  defined  as  x  =  At,,  so  that  the  evaluation  of  the  integral  of  the  product 
of  the  first  term  in  the  summation  of  Eq.  4.21  with  the  Rayleigh  p.d.f.  becomes 


f  OO 

A 

1  T •  f 

2  1 

V  Jo 

— oo 

_  1  / 

1  j 

1  o  \ 

/  n 


8  V  2n.\ 


■'o  Ac 


gfT,  «,  P) 


&(r,  a,  P) 


dr 


(4.28) 


l  I"  -rh/ftr,  «,/?)+! 

Ao 


■dT 


Multiplying  the  second  term  within  the  summation  of  Eq.  4.21  by  the  Rayleigh  p.d.f.  and 
integrating,  yields 
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“■  f))s,{r’  “■  v)™-4*4 

4,tr’  “• <i)' 

r- 

J  /T 


g(x)dxdTfA(A)dA 


(4.29) 


/*T 


4=r  <*,  exp(  -  4pj)g(x)u(jc  -  expj  -  -f—  Ju(,4)rf,4rfx<fl~ 


0  c-oo  */-oo 


v/2tt 


24/ 


where  is  defined  as 


t.-J&nr- «. « 


(4.30) 


Carrying  out  the  integration  with  respect  to  y4  first,  while  realizing  that  w(jc  —  AQu(A) 
is  nonzero  only  for  A  in  the  range  0  <  A  <  yields 


I 


X 

77  A 


~r  exp(  -  A  Pi)  —  exp(  - 
.  2  n 


An 


2  An 


dA 


*0  \  «n  ‘  *  -*  \  £ i 


where  a,  is  defined  as 


ai  =  Pi  + 


1 

24 


(4.31) 


(4.32) 


so  that  Eq.  4.29  becomes 
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A^2n  Jo  Jo  J2n  '  ) 


/  W 

ex  v  £,2 


dxdf 


«  *  A 

-  —  3—  T 

AnJ2n  *’o  •'o 


*o  V 


2aiV'  2?t 


exp  - 


x2  at x7 

2  +  «a 


«? 


(4.33) 


dxdr 


Focusing  on  the  integration  with  respect  to  *  only,  we  have  for  the  first  integral  with 
respect  to  x, 


i 


oo  £ 


-.1 


exp 


jzn 


x 

2 


rf 

”~2^~  exp(  —  ajA2)dAdx 


♦'o 


/  x 
“=■  exp  - 


“'O 
f*  oo 


4at2 


X\a, 


V 


7T  •'O 


f  {'  exp(- 


r2)dr 


dx 


j_  „ 
4)  4(7(-2  v'  2 


exp  - 


.v 


erl 


v/ 


4«pV2 


V 


/ n  /  2  .  -w  4 


and  for  the  second  integral  with  respect  to  x 


/*oo 


‘'O 


2aisj2n 


■  exp 


— =w 


2atylsJ2n 


(4.34) 


(4.35) 


where  y,  is  defined  as 


(4.36) 


2a i 


Thus,  Eq  4.29  becomes 


From  our  previous  definitions  we  can  rewrite  the  expressions  for  a,  and  y,  as 


and 


°‘=P,+TZ 


.  2.4 

Vl  —  I  4 - T~ 


r, «,  0)  +  -^ 

•'0  An 


=  1  + 


~r~gi\ r,  a,  P) 

*'o 


Observe  that  the  terms  in  Eq.  4.37  can  be  expressed  in  terms  of  the  ratio 
signal  energy  to  noise  power  spectral  density  as 


h 

2  -i. 

4/2oa/2 


*(rf  a,  P) 


-h-ffa,  «,  «+l 

i'o 


£<  _ g;(r, «,  /?) 

^(r, 


^druA)dA 

(4.37) 


(4.38) 


(4.39) 


of  average 

(4.40) 

(4.41) 
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so  that  the  bit  error  probability  for  the  third-most  significant  bit  is  obtained  from  the 
sum  of  Eqs.  4.28  and  4.43  -  4.45  .  These  integrals  have  not  been  expressed  in  closed 
form,  however  corresponding  numerical  evaluations  can  and  have  been  carried  out 
without  difficulty.  The  total  BER  under  Rayleigh  fading  conditions  is  finally  given  by 

Pb,  Ray{e)  =  T  pRayi^  incorrect}  +  -J  PRay{ 2B  incorrect} 

,  (4-46) 

+  —  PRay{ 3B  incorrect} 

Appendix  C  presents  a  listing  of  the  FORTRAN  computer  program  used  to  evaluate  the 
BER  under  Rayleigh  fading  conditions  for  the  16-PSK  DBD  receiver.  Data  from  this 
program  is  plotted  and  presented  in  Chapter  5. 

C.  EFFECT  ON  BER  OF  MFSK/QPSK  USING  DBD  RECEIVER 

Recall  that  for  the  MFSK/QPSK  DBD  receiver,  the  BER  is  given  by  Eq.  3.99  as 

This  expression  must  be  rewritten  as  a  function  of  the  amplitude  A  in  order  to  determine 
the  effect  of  Rayleigh  fading  on  receiver  performance.  Substituting  Eq.  4.4  into  Eq.  4.47 
yields 


(4.48) 


Weighting  this  probability  by  the  Rayleigh  p.d.f.  and  integrating  over  the  variable  A,  the 
average  bit  error  probability  in  the  presence  of  signal  fading  is 


Pb,  Ray(e)  =  J°°  c(Aj^yA(A„A 


This  integral  can  be  evaluated  using  the  identity  of  Eq.  4.9  as 


(4.49) 
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^b,  Rayie) 


1 


(4.50) 


D.  EFFECT  ON  BER  OF  MFSK/QPSK  USING  THE  OPTIMUM  RECEIVER 

Evaluation  of  the  (non-fading)  BER  for  MFSK/QPSK  using  the  optimum  receiver 
has  been  carried  out  in  Chapter  3,  Section  C.  Yet  a  different  approach  to  determining 
the  BER  under  Rayeigh  fading  conditions  becomes  more  practical  than  using  modifica¬ 
tions  to  the  approach  previously  used  to  evaluate  the  performance  of  the  DBD  receiver. 
Since  the  BER  given  in  Eq.  3.81  as  a  function  of  s,  the  ratio  of  signal  energy  to  noise 
variance,  namely 


s  = 


(4.51) 


where 

(4.52) 

has  been  previously  evaluated,  a  transformation  of  variables  is  suggested  to  yield 


s  =  g(A)  = 


(4.53) 


where  the  range  of  s  is  0  <  s  S  oo.  Since  the  transformation  is  linear,  the  p.d.f.  of  s  can 
be  written  as 
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where  5^  is  defined  by 


x 


(4.54) 


(4.55) 


The  BER  under  fading  conditions  can  be  evaluated  by  averaging  the  BER  given  by  Eq. 
3.81  over  the  range  of  the  new  random  variable  s  as 

P b,  Ray{e)  ~ 

f  f  C  1  -  +  s)3  N_Vtv)  “fr  exp(  “  TTT  \(s)dyds 

J-ooJ-s  o0  \  2S0  /  (4.56) 

+  (.V  -  1)  f  f  [  1  -  2£>(y)D  N-2[  1  -  Qty  +  S)  -  Q(y  -  s)Dg(V)  —f  exp(  ~  ~rj  V(sKv'<* 

J-ooJ0  S0  \  2 S0  J 


This  integration  also  requires  numerical  techniques,  however  the  program  presented  in 
Appendix  B  can  be  used  as  a  subroutine  for  the  program  listed  in  Appendix  D  which  is 
used  to  evaluate  Eq.  4.56.  Data  from  the  program  of  Appendix  D  is  plotted  and  pre¬ 
sented  in  Chapter  5. 
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V.  RESULTS 

A.  PERFORMANCE  COMPARISONS 

The  symbol  error  rate  (SER)  of  the  receiver  for  M-ary  PSK  signals  is  given  by  Eq. 

1.4,  while  Eq.  1.8  is  used  to  evaluate  the  SER  of  the  receiver  for  M-ary  FSK  signals. 

The  SER  for  MFSK/QPSK  signals  when  using  the  optimum  receiver  of  Figure  12  on 

page  32  is  given  by  Eq.  3.48.  Table  2  lists  the  ratio  of  bit  energv  to  noise  power  spectral 
E 

density  (PSD)  level,  required  to  achieve  symbol  error  rates  of  10~3  and  10-5  for  M- 
A° 

ary  PSK,  M-ary  FSK,  and  MFSK/QPSK  signaling  schemes  with  k  equal  to  3  and  4 

bits'svmbol.  The  table  shows  the  MFSK  QPSK  signaling  scheme  vields  slightlv  better 

E 

performance  than  M-arv  FSK  with  a  nominal  0.2  dB  difference,  in  terms  of  ,  while 

E.  *'o 

outperforming  M-ary  PSK  by  an  average  -rr-  difference  of  5.2  dB  for  k  —  3  and  10.4  dB 

*’0 

for  k  —  4. 


Table  2.  PERFORMANCE  COMPARISONS  AT  KEY  SYMBOL  ERROR 
RATES 


Modulation 

Bits  per 
Symbol 

(dB)  Required  to  Achieve 

iV3  SER  10-‘  SER 

8-PSK 

3 

10.8 

14.4 

16-PSK 

4 

15.5 

IS. 4 

8-FSK 

3 

06.5 

08.6 

16-FSK 

4 

05.6 

07. S 

2 FSK  QPSK 

3 

06.2 

08.5 

4FSK  QPSK 

4 

05.5 

07.6 

The  BER  for  MFSK/QPSK  signaling  using  the  optimum  receiver  shown  in 

Figure  11  on  page  30  is  evaluated  using  Eq.  3.81.  Appendix  B  presents  the  FORTRAN 

computer  program  used  to  compute  this  optimum  receiver's  bit  error  probability  by 

employing  numerical  integration  in  order  to  evaluate  Eq.  3.81.  That  is,  Pb{e}  along  with 

the  union  bound  approximation  to  this  probability  as  given  bv  Eq.  3.82  is  computed  for 
£ 

varying  values  of  -r— .  The  program  is  also  used  to  evaluate  for  the  MFSK  QPSK 


64 


DBD  receiver,  as  given  by  Eq.  3.99.  Data  generated  by  this  program  is  plotted  in 
Figure  16  for  k  =  3,  Figure  17  for  k  =  4,  Figure  18  for  k  =  5,  and  Figure  19  for  k  —  6  . 

In  all  of  the  above  referenced  figures,  numerical  values  for  the  BER  obtained  using 
the  union  bound  approximation  yield  plots  almost  coinciding  with  the  actual  BER  plots 
for  the  optimum  receiver  at  moderate  to  high  -rr-  values.  This  fact  helps  to  confirm  the 
appropriateness  of  the  results  of  the  numerical  integration  used  to  evaluate  the  actual 
receiver  BER.  From  these  plots,  it  becomes  evident  that  the  MFSK/QPSK  DBD  re¬ 
ceiver  suffers  a  performance  penalty  when  compared  to  the  performance  of  the  optimum 

receiver.  Table  3  lists  -rr-  ratios  required  to  achieve  BER's  of  10-3  and  10~6  for  both  the 

*'o 

optimum  and  the  DBD  receiver  structures  for  MFSK/QPSK  modulation.  The  optimum 

receiver  has  an  average  2.5  dB  better  performance  than  the  DBD  receiver  for  k  —  3 

bits  svmbol,  which  increases  to  an  average  value  of  5  dB  for  k  —  4  in  terms  of  -rr-  dif- 

A0 

ference.  For  k  =  6,  the  optimum  receiver  has  an  average  10  dB  advantage  over  the  DBD 
receiver. 


Table  3.  PERFORMANCE  COMPARISONS  OF  MFSK/QPSK  RECEIVERS 


£ 

L'  /  jn\  r» _ • _ i  ^  _  a  _ t  • 


Bits  per 
Svmbol 


y  (dB)  Required  to  Achieve 
10-3  BER  10-6  BER 


05.8 

09.1 

OS.O 

11.8 

05.1 

08.1 

09.8 

13.6 

Modulation' 

Receivers 


2FSK  QPSK 


Opt.  Rcvr 


DBD  Rcvr 


4FSK  QPSK 


Opt.  Rcvr 


DBD  Rcvr 


8FSK  QPSK 


Opt.  Rcvr 


DBD  Rcvr 


8FSK  QPSK 


Opt.  Rcvr 


DBD  Rcvr 


The  reasons  for  this  increased  performance  penalty  can  be  explained  as  follows. 

Eh 

First,  as  can  be  seen  by  examining  the  DBD  receiver  BER  curves  as  a  function  of  -rr-. 

0 

the  performance  of  the  DBD  receiver  degrades  as  k,  the  number  of  bits  per  symbol  is 


04.5 

07.4 

11.8 

15.6 

04.1 

06.8 

14.0 

17.8 
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BIT  ER 
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Figure  17.  4FSK/QPSK  Bit  Error  Performance  (k 


-10  -5  0  5  10  15 


EB/NO  (DB) 


Figure  18.  8FSK/QPSK  Bit  Error  Performance  (k  =  5  bits/symboi) 
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BIT  ERROR  PROBABILITY 


EB/NO  (DB) 


Figure  19.  16FSK/QPSK  Bit  Error  Performance  (k  =  6  bits/symbol) 
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increased.  Recall  from  the  BER  analysis  for  the  DBD  receiver,  that  the  variance  for  the 
summation  output  b„  is  the  sum  of  the  variances  of  the  correlator  outputs  rr  When  the 
number  of  bits  (  and  consequently  the  number  of  correlators  )  is  increased,  there  is  a 
linear  increase  in  the  variance  of  b„  .  This  variance  increase,  which  is  strictly  noise  de¬ 
pendent,  causes  the  bit  error  probability  increase  also.  On  the  other  hand,  examination 
of  the  optimum  receiver  bit  error  curves  reveals  that,  with  increasing  k,  the  optimum 
receiver  exhibits  improved  performance.  This  improvement  reflects  the  tradeoff  that 
exists  between  bit  error  rate  performance  and  signal  bandwidth,  as  increasing  k  results 
in  a  greater  signal  bandwidth.  Note  however  that  this  improved  receiver  performance 
would  not  be  apparent  if  the  bit  error  probability  were  plotted  against  unnormalized 

SNR  instead  of  bit-normalized  SNR,  that  is,  -^r-  Ref.  6  [pp.  172-175]. 

A'o 

B.  BANDWIDTH  EFFICIENCY 

The  bandwidth  efficiency  for  M-ary  PSK,  M-ary  FSK,  and  MFSKQPSK  signaling 
is  given  by  Eqs.  3.105,  3.106  and  3.107,  respectively.  Bandwidth  efficiencies  given  in 
bits  Hz  for  values  of  k  ranging  from  2  to  6  bits/symbol  are  listed  in  Table  4  for  these 
three  signaling  schemes.  The  table  shows  that  M-ary  PSK  is  more  bandwidth  efficient 
than  both  the  M-FSK  and  MFSK  QPSK  schemes.  MFSKQPSK  is  from  three  to  four 
times  as  efficient  as  M-ary  FSK,  with  the  greater  differences  occuring  at  higher  values 
of  k.  MFSK.  QPSK  operates  at  67%  the  efficiency  of  M-ary  PSK  for  k  =  3  bits  symbol, 
and  at  about  40%  the  efficiency  of  M-arv  PSK  for  k  =  4  bits  /symbol. 


Table  4.  BANDWIDTH  EFFICIENCIES  FOR  M-ARY  PSK,  M-ARY  FSK,  AND 
MFSK/QPSK 


Number 
of  Bits 


Number 
of  Sianals 


Bandwidth  Efficiency 

MPSK  MFSK  MFSK  QPSK 


1.00 

0.40 

1.50 

0.33 

2.00 

0.24 

2.50 

0.15 

3.00 

0.09 

C.  EFFECTS  OF  RAYLEIGH  FADING 

The  BER  of  the  8-PSK  DBD  receiver  operating  under  Rayleigh  fading  signal  con- 
ditions  is  given  by  Eq.  4.16  as  a  function  of  average  -rr1  .  This  equation  has  been  used 

•  '  n 


70 


£ 

to  evaluate  receiver  performance  as  a  function  of  average  for  various  values  of  the 
angle  alpha.  The  plots  shown  in  Figure  20  on  page  72  demonstrate  that  the  best  re¬ 
ceiver  performance  occurs  for  a.  =  22.5°.  The  curves  for  a  =  20.0°  and  a  =  25.0°  and  for 
a  =  17.5°  and  a  =  27.5°  are  coincident,  and  therefore  appear  as  one  line  for  each  pair. 

Appendix  C  presents  the  FORTRAN  program  used  to  evaluate  via  numerical  inte¬ 
gration  the  bit  error  probability  under  Rayleigh  fading  signal  conditions  for  the  16-PSK 
DBD  receiver.  Data  generated  by  this  program  is  plotted  in  Figure  21  on  page  73  for 
various  values  of  the  angle  a  while,  maintaining  the  relationship  /?  =  45°  —  a.  This  plot 
shows  the  best  receiver  performance  is  achieved  for  a  =  11.25°  .  Curves  for  a  =  9.75° 
and  ot  =  12.75°  as  well  as  for  a  =  8.25°  and  a.  =  14.25°  are  coincident  and  appear  as  only 
one  line. 


The  receiver's  BER  for  MFSK  QPSK  using  the  DBD  scheme,  is  given  by  Eq.  4.50 
when  Rayleigh  fading  signal  conditions  exist.  Figure  22  on  page  74  illustrates  the  plot 
of  this  BER  for  values  of  k  ranging  from  3  to  5  bits/svmbol.  The  best  performance  is 
obtained  for  k  equal  to  3,  and  performance  degrades  with  increasing  values  of  k.  The 
performance  curve  has  a  near  linear  characteristic  for  values  of  average  ~~  exceeding 
about  5  dB  as  opposed  to  the  "waterfall"  characteristic  of  the  BER  curves  when  Rayleigh 
fading  conditions  were  not  present. 

The  BER  for  MFSK  QPSK  signaling  using  the  optimum  receiver  under  Rayleigh 
fading  signal  conditions  is  given  by  Eq.  4.56  .  Appendix  D  presents  the  listing  of  the 
FORTRAN  program  used  to  evaluate  this  BER,  which  involves  double  numerical  inte¬ 
gration.  Data  generated  by  this  program  is  plotted  in  Figure  23  on  page  75  for  values 
of  k  ranging  from  3  to  5  bits  ' symbol.  The  curves  show  a  somewhat  surprising  result  in 
that  the  best  performance  is  obtained  when  k  is  equal  to  3,  and  performance  degrades 
with  increasing  values  of  k.  This  result  is  opposite  to  that  encountered  when  Rayleigh 
fading  conditions  were  not  present,  in  which  increasing  k  produced  improved  receiver 


performance.  The  curves  of  Figure  23  also  exhibit  nearly  linear  characteristics  for  val¬ 
ues  of  average  ~~  greater  than  approximately  0  dB. 

Table  5  on  page  76  presents  a  summary  of  the  average  values  required  to 
achieve  bit  error  rates  of  KM  and  KM  for  all  of  the  receivers  discussed  above,  under  both 


non-fading  and  Rayleigh  fading  signal  conditions.  For  the  MFSK/QPSK  receivers,  in 

order  to  simplifiy  the  comparisons,  only  the  results  for  k  =  3  are  presented.  For  the 

MFSK  QPSK  receivers,  the  performance  degradation  is  about  the  same  for  both  the 

optimum  and  the  DBD  receivers.  In  order  to  achieve  a  BER  of  1CM  both  receivers  suffer 

a  nominal  3  dB  performance  loss  in  terms  of  average  under  Rayleigh  fading  condi- 

^'0 
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BIT  ERROR  PROBAB 


Figure  20.  8-PSK  DBD  Receiver  Bit  Error  Performance  in  Rayleigh  Fading 
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BIT  ERROR  PROBABILITY 
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BIT  ERROR  PROBABILITY 


AVERAGE  EB/NO  (DB) 

Figure  22.  MFSK/QPSK  DBD  Rcvr.  Bit  Error  Performance  in  Rayleigh  Fading 


74 


LEGEND 
K  =  3 


tions,  and  performance  is  correspondingly  degraded  by  about  10  dB  at  a  BER  of  10~2  . 
For  the  M-ary  PSK  DBD  receivers,  performance  is  degraded  by  about  4  dB  and  by  12 
dB  at  a  BER  of  10-'  and  10~2,  respectively. 


Table  5.  EFFECT  OF  RAYLEIGH  FADING  USING  VARIOUS  RECEIVERS 


Modulation/ 

Receiver 

Rayleigh 

Fading 

Eb 

— (dB)  Required  to  Achieve 
l°0-»  BER  10-2  BER 

2FSK  QPSK 

Optimum 

Receiver 

No 

-0.6 

03.7 

Yes 

02.9 

14.0 

2FSK  QPSK 

DBD 

No 

00.4 

05.6 

Receiver 

Yes 

03.6 

15.1 

S-PSK 

DBD 

No 

00.1 

03.9 

Receiver 

Yes 

04.4 

16.0 

16-PSK 

DBD 

No 

03.6 

07.4 

Receiver 

Yes 

07.2 

19.4 
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VI.  CONCLUSIONS 


The  results  presented  show  that  MFSK/QPSK  modulation  is  superior  to  M-ary  FSK 
in  aspects  of  error  performance  and  bandwidth  efficiency.  While  MFSK/QPSK  has  only 
slightly  better  error  performance  over  M-ary  FSK,  MFSK/QPSK  requires  only  one-third 
to  one-fourth  the  bandwidth  used  by  M-ary  FSK  signaling.  Also,  since  the 
MFSK/QPSK  receiver  requires  the  use  of  only  one  fourth  the  number  of  synchronized 
carriers  and  integrate  and  dump  circuits  required  by  the  M-ary  FSK  receiver,  it  can  be 
argued  that  the  MFSK/QPSK  receiver  is  less  complex  than  the  M-ary  FSK  receiver. 
The  results  also  show  that  MFSK/QPSK  is  comparable  to  M-ary  PSK  in  the  sense  that 
it  offers  improved  error  performance  over  M-ary  PSK,  but  does  not  reduce  bandwidth 
efficiency  as  much  as  orthogonal  signaling. 

In  comparing  the  different  types  of  MFSK/QPSK  receivers,  the  DBD  receiver  offers 
the  ability  to  decode  bits  individually  without  symbol-to-bit  mapping  circuitry;  however, 
this  is  accomplished  at  the  expense  of  error  performance.  This  error  performance  pen¬ 
alty  becomes  prohibitive  for  larger  values  of  k.  Since  the  advantage  of  receiver  simplicity 
is  no  longer  as  attractive  for  larger  k.  the  DBD  receiver  should  be  considered  as  an  al¬ 
ternative  to  the  optimum  receiver  only  for  smaller  values  of  k. 

Analysis  of  the  effect  of  Rayleigh  fading  in  M-ary  FSK  DBD  receivers  and  in 
MFSK  QPSK  receivers  has  generally  demonstrated  severe  performance  degradation  in 
the  sense  of  producing  large  increases  in  the  bit  energy  to  noise  PSD  ratios  required  to 
maintain  only  modest  BER  performance.  The  penalties  increase  rapidly  for  lower  bit 
error  probabilities.  Unless  fading  effects  can  be  overcome,  all  of  the  receivers  studied 
become  impractical  in  the  presence  of  fading.  The  MFSK  QPSK  receivers  show  a  slight 
reduction  in  performance  degradation  over  the  M-ary  PSK  DBD  receivers.  For  the 
M-ary  PSK  receivers,  the  same  value  of  the  angle  adjustment  a.  which  produced  the 
minimum  BER  for  a  non-fading  channel  also  produced  the  best  overall  performance 
under  Rayleigh  fading  conditions. 
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APPENDIX  A.  INTEGRAL  IDENTITY 

This  appendix  presents  the  proof  for  the  integral  identity  used  in  the  text  as  Eq.  4. 13. 


f  °°  Q{<*A)Q(pAyx{A)dA  = 


1 

4 


P  tan 


-1 


'p2  +  J T 

Al 


a 


2N*2+4 


(5.1) 


We  integrate  by  parts,  i.e., 


udv  =  uv 


with  the  following  substitutions 


u  =  Q(<xA)Q(PA) 


and 


dv  =fx(A)dA 


Differentiating  Eq.  B.3  with  respect  to  A ,  yields 

~  =  Q(*A)Q(PA)  4-  QizA)Q(PA) 


(B.2) 


(B-3) 


(BA) 


(5.5) 


where 
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(*  oo 


_d_ 

dA 


7==-exp(  -V  )dy 

;2  n 


(5.6) 


i2n 


exp 


2  a  2 
a  /j 


where  the  last  equality  is  obtained  using  Liebnitz  rule  to  differentiate  the  integral.  Sim¬ 
ilar  results  hold  for  Q{pA),  thus  yielding 


du  —  — 


Q(aA)  —jLr  expf  -  -  'j  +  Q(pA) -jL=r  exp^  - 

V  L  J  ^2n  \  <- 


dA  ( B.l ) 


Now,  integrating  Eq.  B.4  to  obtain  v  results  in 


J  Ao 


A  exp  (  -  )dA 


2  At, 


—\~14 


Solving  for  the  first  term  on  the  right  hand  side  of  Eq.  B.2  yields 


i 

Wn  = 


-Q(*A)Q(pA)  exp  - 


A‘ 


=  Q(0) 


X 

4 


For  the  second  term  on  the  right  hand  side  of  Eq.  B.2  ,  we  have 


(5.8) 


(5.9) 


J  vdu  =  j  ^(-“7 

Jo  Jo  \  2  A0 


Q(aA)-=-exp(-~—  )  +  Q(PA)  -  5=r  expl  - 


'  2  JT 


J2n 


2  x2 

a  /! 


(5.10) 


/»oor 


CM)1 


-TT 


exp(  -  z2A2)  +  Q{pA)-  —  exp(  -  z2/!2) 


•'o 
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where  z,  and  z2  are  defined  as 


and 


(5.11) 


(5.12) 


The  form  of  these  integrals  is  similar  to  the  identity  given  by  [Ref.  9:  pp.  1-20]  as 


f—  'xp - 1^7  ln[  -frf  ] 

where  a  and  b  may  be  complex,  and 


_1 

Real(a)  "1 

tan 

Imag(a) 

<i 


Using  this  identity  and  the  equality 


(?(•*)  =  *2"  Erfcf  — = 


produces  an  equivalent  form,  namely 

l 

Using  Eq.  B.  16  to  solve  Eq.  B.10  yields 


Q(ax)  exp(b2x2)dx  =  '' 


In 


4/bN:  n 


q  +  b^’2 
a  -  bj2 


(5.13) 


(5.14) 


(5.15) 


(5.16) 


j4nzx!2 


«  +jz  iV2 

«  -jz iv'2 


+  • 


jAnZyjl 


-=-  In 


P  ±j± 2V'2 
P  ~J:2\  2 


(5.17) 


Focusing  on  the  argument  of  the  first  log  function  and  converting  to  complex  exponen¬ 
tial  notation  produces 


a  +jzjl_ 
a  -jzxyj2 


sj  a2  4-  2 z\  exp 


^ !"'2 


v/a2  +  2212  exp 


i«»'<  ~^2 


=  exp 


=  exp 


tan  H  -1—-  )  -  tan  1 


'y2tan-.fiiA 


(5.18) 


where  the  last  equality  holds  by  restricting  both  a  and  z,  to  be  positive,  which  is  indeed 
the  case.  Using  this  result  in  Eq.  B.17  yields 


p  tan-1 


2n 


+ 


a  tan 


(5.19) 


Finally,  substituting  Eqs.  B.9  and  B.19  into  Eq.  B.2  yields  the  right  hand  side  of  Eq.  B.l, 
thus  proving  the  identity. 
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APPENDIX  B.  MFSK/QPSK  BER  COMPUTER  PROGRAM 


* 

*  PROGRAM  TO  COMPUTE  THE  BIT  ERROR  RATES  FOR  THE  BI -ORTHOGONAL 

*  SIGNAL  SET  (BOSS)  FOR  THE  OPTIMUM  RCVR  AND  FOR  THE  MFSK/QPSK 

*  DIRECT  BIT  DETECTION  (DBD)  RCVR. 

* 

* 

************************************************************************ 


PROGRAM  PRBTER 


IMPLICIT  CHARACTER*!  (A-Z) 


DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
INTEGER 
INTEGER 


EBNODB,  INIT, FINAL, DLTADB 
EBDVNO , SREDSI 

PRBTCR , PRBERR , INTGL1 , INTGL2 , INTGL3 
DBDPBE .UBAPBE 

DCADRE , F 1 , A 1 , B 1 , AERR , RERR , ERROR 
F2,A2,B2,F3 
PIE,SQR2PI ,SQR8PI,0NE 
I , K , N , M , F , IER , NUMPTS 
LEVEL, LEVOLD 


EXTERNAL  F1,F2,F3 

COMMON  SREDSI ,N 


C 


INITIALIZE  CONSTANTS 


PIE 

= 

3. 141592654D0 

SQR2PI 

= 

DSQRT( 2*PIE) 

SQR8PI 

= 

DSQRTC  8*PIE) 

K 

= 

6 

M 

= 

2**K 

N 

= 

2**(K-1) 

F 

= 

2**(K-2) 

A2 

= 

00. 0D0 

Bl 

= 

20. 0D0 

B2 

= 

20. 0D0 

ONE 

= 

01.  0D0 

NUMPTS 

= 

101 

INIT 

= 

-10. 0D0 

FINAL 

= 

8.  0D0 

DLTADB 

= 

(FINAL  -  INIT) /(NUMPTS- 1) 

AERR 

= 

0.  0 

RERR 

= 

. 0000001 

C 

LEVEL  = 


SUPPRESS  ALL  ERROR  MESSAGES. 
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CALL  UERSET( LEVEL, LE VOID) 
CALL  XUFLOW(O) 


COMPUTE  PROBABILITY  OF  BIT  ERROR 
FOR  VALUES  OF  EB/NO  IN  THE  RANGE  OF: 
INIT  <  EBNODB  <  FINAL 

WRITE(8,100)  K,F 
WRITE(8,*)'  ' 

WRITE(8,*)'  OPTIMUM  RCVR  DBD  RCVR 

&  UNION  BOUND’ 

WRITE(8,*)'  BIT  ERROR  BIT  ERROR 

&  APPROXIMATION' 

WRITE(8 ,*) ' INDEX  EB/NO  (DB)  PROBABILITY  PROBABILITY 

&  BIT  ERROR’ 

EBNODB  =  INIT 
DO  10  I=1,NUMPTS 

EBDVNO  =  10. 0D0**( EBNODB/ 10) 

SREDSI  =  DSQRT( 2*£BDVN0*K) 

A1  =  -SREDSI 

INTGL1  =  DCADRE (FI, A1,B1,AERR,RERR, ERROR, IER)  /  SQR2PI 
INTGL2  =  DCADRE(F2,A2,B2,AERR,RERR, ERROR, IER)  /  SQR8PI 
INTGL3  =  DCADRE(F3,A1,B2,AERR,RERR, ERROR, IER)  /  SQR8PI 
IF( INTGL3.  GT.  0.  ODO)  THEN 

PRBTCR  =  (N-l)*( INTGL2+INTGL3)  +  INTGL1 

ELSE 

PRBTCR  =  (N-1)*INTGL2  +  (N-1)*INTGL3  +  INTGL1 

END  IF 

PRBERR  =  ONE  -  PRBTCR 

DBDPBE  =  DERFC(  DSQRT( 2*K*EBDVN0/M)  )  /  2 
UBAPBE  =  (  (M  -  2)*DERFC(  DSQRT(K*EBDVN0/2)  ) 

&  +  DERFC(  DSQRT( K*EBDVNO )  )  )  /U 

WRITE( 8,101)  I .EBNODB , PRBERR .DBDPBE , UBAPBE 
EBNODB  =  EBNODB  +  DLTADB 
10  CONTINUE 


100  F0RMAT(2X,' NUMBER  OF  BITS: ', 2X, I 1 , 6X, ' NUMBER  OF  FREQUENCIES: ', 

5c  2X.I1) 

101  FORMAT( 2X, 13 , 3X,F10.  5 , 3X,E14.  5E2,3X,E14. 5E2,3X,E14. 5E2) 

STOP 

END 


DOUBLE  PRECISION  FUNCTION  F1(Y) 

DOUBLE  PRECISION  Y, SREDSI , SQRT2 
INTEGER  N 

COMMON  SREDSI ,N 

C  INITIALIZE  CONSTANT 

SQRT2  =  DSQRT(2. ODO) 
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C  COMPUTE  F1(Y,SREDSI,N) 

FI  =  DERF(  (Y+SREDSI)/SQRT2  )**(N-1)  *  DEXP( -(Y**2)/2) 

RETURN 

END 


DOUBLE  PRECISION  FUNCTION  F2(Y) 

DOUBLE  PRECISION  Y , FACTR 1 , FACTR2 , SREDS I 
DOUBLE  PRECISION  SQRT2 
INTEGER  N 


COMMON  SREDS I, N 

C  INITIALIZE  CONSTANTS 

SQRT2  =  DSQRTC2.  ODO) 

C  COMPUTE  F2(Y,SREDSI,N) 

FACTR1  =  DERF(  Y/SQRT2  )**(N-2) 

FACTR2  =  DERF(  (Y+SREDSI)/SQRT2  ) 

F2  =  FACTRl*FACTR2*DEXP(-(Y**2)/2) 

RETURN 

END 


************************************************************************ 
DOUBLE  PRECISION  FUNCTION  F3(Y) 

DOUBLE  PRECISION  Y,FACTR1 ,FACTR2 ,SREDSI 
DOUBLE  PRECISION  SQRT2 
INTEGER  N 


COMMON  SREDS I ,N 


C 


SQRT2  =  DSQRT(2. ODO) 


INITIALIZE  CONSTANTS 


C 


COMPUTE  F3(Y,SREDSI,N) 


FACTR 1  =  DERF(  (Y+SREDSI)/SQRT2  )**(N-2) 

FACTR2  =  DERF(  Y/SQRT2  ) 

F3  =  F ACTR 1*FACTR2*DEXP( - ( Y+SREDS I ) **2 / 2  ) 


RETURN 

END 
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APPENDIX  C.  16-PSK  RAYLEIGH  FADING  BER  COMPUTER 

PROGRAM 


PROGRAM  DBDRAY 


IMPLICIT  CHARACTER* 1  (A-Z) 

DOUBLE  PRECISION  EBNODB , INIT, FINAL, DLTADB ,PBERAY 
DOUBLE  PRECISION  EBDVNO,AVESNR, ALPHA, AL$DEG, BETA 
DOUBLE  PRECISION  P1BERY ,V1$FTN(4) ,ARGMNT 
DOUBLE  PRECISION  P2BERY,V2$FTN(4) ,ARG$1 ,ARG$2 
DOUBLE  PRECISION  P3BERY, INTGRL( 8) ,TERM$ 1 ,TERM$2 ,TERM$3 
DOUBLE  PRECISION  FI , A,B ,ERRABS .ERRREL, INTGL1 ,ERREST 
DOUBLE  PRECISION  F2 ,F3 , INTGL2, INTGL3 
DOUBLE  PRECISION  PI ,ONE,SUM$l ,SUM$2 
INTEGER  I , INDEX , NUMPTS , IRULE 

EXTERNAL  F1,F2,F3 

COMMON  AVESNR, ALPHA, BETA, INDEX 


C  SUPPRESS  IMSL10  ERROR  MESSAGES. 

CALL  ERSET( 3 ,0 , -1) 

C  INITIALIZE  CONSTANTS 

PI  =  3. 141592654D0 

ONE  =1. 0D0 
AL$DEG  =  14. 25D0 

ALPHA  =  AL$DEG  *  PI/ 180 
BETA  =  PI/4  -  ALPHA 

V1$FTN( 1)  =  DCOS( ALPHA) 

V1$FTN( 2)  =  DCOS(BETA) 

V1$FTN( 3)  =  DSIN( ALPHA) 

V1$FTN(4)  =  DSIN(BETA) 

V2$F'iN(l)  =  DCOS( ALPHA)  -  DSIN( ALPHA) 

V2§FTN( 2)  =  DCOS( ALPHA)  +  DSIN( ALPHA) 

V2$FTN( 3)  =  DCOS( BETA)  -  DSIN( BETA) 
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V2$FTN(4)  =  DCOS(BETA)  +  DSIN(BETA) 


NUMPTS 

=5 

101 

INIT 

= 

-10. ODO 

FINAL 

= 

20. ODO 

DLTADB 

= 

(FINAL  -  INIT)/( NUMPTS- 1) 

A 

s 

0.  ODO 

B 

= 

PI/8 

ERRABS 

= 

0.  ODO 

ERRREL 

= 

0. 0001 

IRULE 

=5 

2 

COMPUTE  PROBABILITY  OF  BIT  ERROR 
FOR  VALUES  OF  EBNODB  IN  THE  RANGE  OF: 
INIT  <  EBNODB  <  FINAL 


WRITE( 8,*)' EFFECT  OF  RAYLEIGH  FADING  ON  16-PSK,  DBD  RCVR* 
WRITE( 8 , 100)  AL$DEG 
WRITE( 8 ,*) '  ' 

WRITEC8,*)'  DBD  RCVR' 

WRITEC8,*)'  BIT  ERROR' 

WRITE(8,*)' INDEX  EB/NO  (DB)  PROBABILITY' 

EBNODB  =  INIT 


C  BEGIN  MAIN  LOOP. 

DO  70  1=1 , NUMPTS 
EBDVNO  =  10.0D0**(EBN0DB/10) 

AVESNR  =  2  *  EBDVNO 

C  CALCULATE  PROBABILITY  OF  IB  ERROR 

SUM$1  =  0. ODO 

DO  10  INDEX=1 ,4 

ARGMNT  =  EBDVNO  /  (4  *  EBDVNO  *  V1$FTN( INDEX)**2  +  ONE) 
SUM$1  *  SUM$ 1  +  V1$FTN( INDEX)  *  DSQRT( ARGMNT) 

10  CONTINUE 

P1BERY  =  (2. ODO  -  SUM$l)/4 

C  CALCULATE  PROBABILITY  OF  2B  ERROR 

SUM$1  =  0.  ODO 

DO  20  INDEX=1 ,4 

ARG$1  =  EBDVNO  /  (EBDVNO  *  V2$FTN( INDEX )**2  +  ONE/2) 

SUM$1  =  SUM$1  +  V2$FTN( INDEX)  *  DSQRT(ARG$1) 

20  CONTINUE 


SUM$2  =  0. ODO 

DO  30  INDEX=1 ,4 

ARG$1  =  DSQRT(  V2$FTN( INDEX)**2  +  ONE  /  (2*EBDVN0)  ) 
ARG$2  =  ARG$ 1  /  V2$FTN(  INDEX  -  ( -1)**( INDEX)  ) 

SUM$2  =  SUM$2  +  V2$FTN( INDEX)  *  DATAN(ARG$2)  /  ARG$1 
30  CONTINUE 

P2BERY  =  ONE/2  -  SUM$l/4  +  SUM$2/(2*PI) 

C  CALCULATE  PROBABILITY  OF  3B  ERROR 
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DO  40  INDEX=1 , 8 

CALL  DQDAG  (FI ,A,B .ERRABS .ERRREL, IRULE , INTGL1 .ERREST) 
TERM$1  =  DSQRT( AVESNR ) * I NTGL 1/4 

CALL  DQDAG  ( F2 , A , B , ERRAB S , ERRREL .IRULE, INTGL2 , ERREST ) 
TERM$2  =  INTGL2 

CALL  DQDAG  ( F3, A, B, ERRABS , ERRREL, IRULE, INTGL3, ERREST) 
TERM$3  =  DSQRT( AVESNR )*INTGL3 

INTGRL( INDEX)  =  TERM$1  -  (TERM? 2  +  TERM$3)/(4*PI) 

40  CONTINUE 

SUM$1  =  0.  0D0 

DO  50  INDEX=5 ,8 
SUM$1  =  SUM? 1  +  INTGRL( INDEX) 

50  CONTINUE 

SUM?2  =  O.ODO 

DO  60  INDEX® 1,4 
SUM $2  =  SUM $2  +  INTGRLC INDEX) 

60  CONTINUE 

P3BERY  =  ONE/2  +  SUM?1  -  SUM$2 
PBERAY  =  P1BERY/2  +  P2BERY/4  +  P3BERY/4 
WRITE(8, 101)  I, EBNODB, PBERAY 
EBNODB  =  EBNODB  +  DLTADB 
70  CONTINUE 

C  END  MAIN  LOOP 

100  F0RMAT(2X, 'ALPHA  =' , IX ,F5. 2 , IX, ' DEGREES.  ' ) 

101  F0RMAT(2X,I3,3X,F10.  5,3X,E14.  5E2,3X,E14. 5E2,3X,E14. 5E2) 

STOP 

END 


DOUBLE  PRECISION  FUNCTION  FI (GAMMA) 


* 


IMPLICIT  CHARACTER* 1  (A-Z) 

DOUBLE  PRECISION  GAMMA .AVESNR, ALPHA, BETA 
DOUBLE  PRECISION  TERM? 1 , FINDEX , G INDEX 
INTEGER  INDEX 

COMMON  AVESNR , ALPHA , B  ETA , I NDEX 


COMPUTE  Fl( INDEX , GAMMA , AVESNR , ALPHA , BETA) 

GO  TO  (10,20,30,40,50,60,70,80)  INDEX 

10  FINDEX  =  DSIN( GAMMA  -  ALPHA) 

GINDEX  =  DCOS( GAMMA  -  ALPHA) 

GO  TO  100 


87 


o  o 


20  FINDEX  =  DCOS( GAMMA  -  ALPHA) 

G INDEX  =  DSIN( GAMMA  -  ALPHA) 

GO  TO  100 

30  FINDEX  =  DCOS( GAMMA  +  ALPHA) 

G INDEX  =  DSIN( GAMMA  +  ALPHA) 

GO  TO  100 

40  FINDEX  =  DSIN( GAMMA  +  ALPHA) 

GINDEX  =  DCOS( GAMMA  +  ALPHA) 

GO  TO  100 

50  FINDEX  =  DSIN( GAMMA  -  BETA) 

GINDEX  =  DCOS( GAMMA  -  BETA) 

GO  TO  100 

60  FINDEX  =  DCOS( GAMMA  -  BETA) 

GINDEX  =  DSIN( GAMMA  -  BETA) 

GO  TO  100 

70  FINDEX  =  DCOS( GAMMA  +  BETA) 

GINDEX  =  DSIN( GAMMA  +  BETA) 

GO  TO  100 

80  FINDEX  =  DSIN( GAMMA  +  BETA) 

GINDEX  =  DCOS( GAMMA  +  BETA) 

100  TERM? 1  =  AVESNR  *  FINDEX**2  +  1. 0D0/2 
FI  =  GINDEX  /  TERM?1  /  DSQRT(TERM$1) 

RETURN 

END 


j.  -V 


^  yrVf  V*  "tC  VrVf  Vc  Vc  #V 1 


DOUBLE  PRECISION  FUNCTION  F2( GAMMA) 


IMPLICIT  CHARACTER* 1  (A-Z) 

DOUBLE  PRECISION  GAMMA .AVESNR, ALPHA .BETA 
DOUBLE  PRECISION  TERM? 1 ,TERM$2 .FINDEX, GINDEX 
INTEGER  INDEX 

COMMON  AVESNR, ALPHA, BETA, INDEX 


COMPUTE  F2 ( INDEX , GAMMA , AVESNR , ALPHA , BETA) 


GO  TO  (10,20,30,40,50,60,70,80)  INDEX 

10  FINDEX  =  DSIN( GAMMA  -  ALPHA) 

GINDEX  =  DCOS( GAMMA  -  ALPHA) 

GO  TO  100 


20  FINDEX  =  DCOS( GAMMA  -  ALPHA) 
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GINDEX  =  DSIN( GAMMA  -  ALPHA) 

GO  TO  100 

30  FINDEX  =  DCOS( GAMMA  +  ALPHA) 

GINDEX  =  DSIN( GAMMA  +  ALPHA) 

GO  TO  100 

40  FINDEX  =  DSIN( GAMMA  +  ALPHA) 

GINDEX  =  DCOS( GAMMA  +  ALPHA) 

GO  TO  100 

50  FINDEX  =  DSIN( GAMMA  -  BETA) 

GINDEX  =  DCOS( GAMMA  -  BETA) 

GO  TO  100 

60  FINDEX  =  DCOS( GAMMA  -  BETA) 

GINDEX  =  DSIN( GAMMA  -  BETA) 

GO  TO  100 

70  FINDEX  =  DCOS( GAMMA  +  BETA) 

GINDEX  =  DSIN( GAMMA  +  BETA) 

GO  TO  100 

80  FINDEX  =  DSIN( GAMMA  +  BETA) 

GINDEX  =  DCOS( GAMMA  +  BETA) 

100  TERM? 1  =  AVESNR  *  FINDEX**2  +  1.  0D0/2 

TERM? 2  =  AVESNR  *  GINDEX**2 
F2  =  TERM? 2  /  (TERM? 1*TERM?2  +  TERM? 1**2) 

RETURN 

END 


DOUBLE  PRECISION  FUNCTION  F3( GAMMA) 


'  ■'« 


IMPLICIT  CHARACTER* 1  (A-Z) 

DOUBLE  PRECISION  GAMMA .AVESNR, ALPHA .BETA 

DOUBLE  PRECISION  TERM? 1 ,TERM?2 .FACTRl .FACTR2 .FINDEX, GINDEX 

INTEGER  INDEX 

COMMON  AVESNR, ALPHA, BETA, INDEX 


COMPUTE  F3( INDEX , GAMMA , AVESNR , ALPHA , BETA ) 

GO  TO  (10,20,30,40,50,60,70,80)  INDEX 

10  FINDEX  =  DSIN( GAMMA  -  ALPHA) 

GINDEX  =  DCOS( GAMMA  -  ALPHA) 

GO  TO  100 

20  FINDEX  =  DCOS( GAMMA  -  ALPHA) 
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G INDEX  =  DSIN( GAMMA  -  ALPHA) 

GO  TO  100 

30  FINDEX  =  DCOS( GAMMA  +  ALPHA) 

GINDEX  =  DSIN( GAMMA  +  ALPHA) 

GO  TO  100 

40  FINDEX  =  DSIN( GAMMA  +  ALPHA) 

GINDEX  =  DCOS( GAMMA  +  ALPHA) 

GO  TO  100 

50  FINDEX  =  DSIN( GAMMA  -  BETA) 

GINDEX  =  DCOS( GAMMA  -  BETA) 

GO  TO  100 

60  FINDEX  =  DCOS( GAMMA  -  BETA) 

GINDEX  =  DSIN( GAMMA  -  BETA) 

GO  TO  100 

70  FINDEX  =  DCOS( GAMMA  +  BETA) 

GINDEX  =  DSIN( GAMMA  +  BETA) 

GO  TO  100 

80  FINDEX  =  DSIN( GAMMA  +  BETA) 

GINDEX  =  DCOS( GAMMA  +  BETA) 

100  TERMS  1  =  AVESNR  *  FINDEX**2  +  1. ODO 

TERM$ 2  =  AVESNR  *  GINDEX**2 
FACTR1  *  GINDEX  /  TERM$1  /  DSQRT(TERM$1) 
FACTR2  =  DATAN(  DSQRT(TERM$2/TERM$1)  ) 

F3  =  FACTR1  *  FACTR2 


RETURN 

END 
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APPENDIX  D.  MFSK/QPSK  RAYLEIGH  FADING  BER  COMPUTER 

PROGRAM 


* 

* 

*  PROGRAM  TO  COMPUTE  THE  BIT  ERROR  RATES  FOR  THE  B I -ORTHOGONAL 

*  SIGNAL  SET  (BOSS)  FOR  THE  OPTIMUM  RCVR  WHEN  SUBJECTED  TO  RAYLEIGH 

*  FADING. 

* 

* 


PROGRAM  PBTERA 


IMPLICIT  CHARACTER* 1  (A-Z) 

DOUBLE  PRECISION  EBNODB , INIT, FINAL, DLTADB 
DOUBLE  PRECISION  EBDVNO ,S0SQRD,S 
DOUBLE  PRECISION  PBCRAY , PBERAY 

DOUBLE  PRECISION  DCADR2 , INTGND,LIL,UIL,AERR,RERR, ERROR 
DOUBLE  PRECISION  ONE 
INTEGER  I ,K,N,M,F,IER,NUMPTS 

INTEGER  LEVEL ,LEVOLD 

EXTERNAL  INTGND 

COMMON  SOSQRD.S.N 


C  INITIALIZE  CONSTANTS 

K  =5 

M  =  2**K 

N  =  2**(K-1) 

F  =  2**(K-2) 

LIL  =  00. 0D0 

ONE  =  01. 0D0 

NUMPTS  =  101 

INIT  =  -10. 0D0 
FINAL  *  20. 0D0 

DLTADB  =  (FINAL  -  INIT)/(NUMPTS-1) 

AERR  =0.0 
RERR  =  .0001 

C  SUPPRESS  ALL  ERROR  MESSAGES. 

LEVEL  =  0 

CALL  UERSET( LEVEL, LEVOLD) 

C  CALL  XUFLOW ( 0 ) 
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COMPUTE  PROBABILITY  OF  BIT  ERROR 
FOR  VALUES  OF  EB/NO  IN  THE  RANGE  OF: 
INIT  <  EBNODB  <  FINAL 

WRITE( 8 ,*)’ EFFECT  OF  RAYLEIGH  FADING  ON  BOSS’ 

WRITE(8,*)'  ’ 

WRITE( 8 , 100)  K,F 
WRITE(  8  ,*) '  ’ 

VRITE( 8 ,*) '  OPTIMUM  RCVR’ 

WRITE(8,*)'  BIT  ERROR  ' 

WRITE( 8  ,*) ' INDEX  EB/NO  (DB)  PROBABILITY  ' 

EBNODB  =  INIT 
DO  10  1=1 ,NUMPTS 

EBDVNO  =  10.  0D0**(EBN0DB/10) 

SOSQRD  =  EBDVNO  *  K 

UIL  =  MIN( 18.  5D0,10*DSQRT( SOSQRD)  ) 

PBERAY  =  DCADR2(INTGND, LIL.UIL.AERR.RERR, ERROR, IER)  /  SOSQRD 
WRITE( 8 , 101)  I, EBNODB, PBERAY 
EBNODB  =  EBNODB  +  DLTADB 
10  CONTINUE 

100  F0RMAT(2X,' NUMBER  OF  BITS:  ' , 2X, II, 6X, 'NUMBER  OF  FREQUENCIES:  ' , 

&  2X , I 1 ) 

101  FORMAT( 2X, 13 ,3X,F10.  5 , 3X,E14.  5E2,3X,E14.  5E2,3X,E14.  5E2) 

STOP 

END 


****************iV*********Vr****Vr**************************************** 

DOUBLE  PRECISION  FUNCTION  INTGND(SREDSI) 


IMPLICIT  CHARACTER*!  (A-Z) 


DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
INTEGER 


EBDVNO , SREDSI .SOSQRD , S .THRHLD 
PRBTCR , PRBTER , ORTERM , APTERM 
INTGL1 , INTGL2 , INTGL3 
DC ADRE , F 1 , A 1 , B 1 , AERR , RERR , ERROR 
F2,A2,B2,F3 
PIE , SQR2PI , SQR8PI , ONE 
N,  IER 


EXTERNAL 

COMMON 


F1,F2,F3 
SOSQRD, S,N 


THRHLD  =  6.  2D0 
PIE  =  3. 141592654D0 
ONE  =  l.OODO 
SQR2PI  =  DSQRT( 2*PIE) 
SQR8PI  =  DSQRT( 8*PIE) 
AERR  =0.0 
RERR  =  .0000001 


INITIALIZE  CONSTANTS 


o  o 


A2 

=  00. 0D0 

B1 

=  20. ODO 

B2 

=  20. ODO 

COMPUTE  PROBABILITY  OF  BIT  ERROR  AS  A 
FUNCTION  OF  SREDSI  AND  N 

IF(SREDSI.GT.THRHLD)  THEN 
C  USE  UNION  BOUND  APPROXIMATION 

ORTERM  =  (2*N  -  2)*DERFC(  SREDSI/2  ) 

APTERM  ■  DERFCC  SREDSI  /  DSQRT(2.0D0)  ) 

PRBTER  =  (  ORTERM  +  APTERM  )  /  4 
ELSE 

C  USE  EXACT  CALCULATION 

S  =  SREDSI 

A1  =  -SREDSI 

INTGL1  =  DCADRE ( F 1 , A1 , B 1 , AERR , RERR , ERROR ,  IER )  /  SQR2PI 
INTGL2  =  DCADRE(F2 ,A2 ,B2 , AERR, RERR, ERROR, IER)  /  SQR8PI 
INTGL3  =  DCADRE(F3,A1,B2, AERR, RERR, ERROR, IER)  /  SQR8PI 
IF( INTGL3. GT. 0.  0D0)  THEN 

PRBTCR  ■  INTGL1  +  (N-l)*( INTGL2+INTGL3) 

ELSE 

PRBTCR  =  INTGL1  +  (N-1)*INTGL2  +  (N-1)*INTGL3 

ENDIF 

PRBTER  =  ONE  -  PRBTCR 


ENDIF 


C  COMPUTE  INTGNDC SREDSI , SOSQRD ,N) . 

INTGND  =  PRBTER  *  SREDSI  *  DEXP(  “(SREDSI**2)  /  (2  *  SOSQRD)  ) 

RETURN 

END 


************************************************************************ 
DOUBLE  PRECISION  FUNCTION  F1(Y) 

DOUBLE  PRECISION  Y, SOSQRD, S.SQRT2 
INTEGER  N 

COMMON  SOSQRD, S,N 

C  INITIALIZE  CONSTANT 

SQRT2  =  DSQRT(2.0D0) 

C  COMPUTE  F1(Y,S,N) 

FI  =  DERF(  ( Y+S ) / SQRT2  )**(N-1)  *  DEXP( -(Y**2)/2) 

RETURN 

END 
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DOUBLE  PRECISION  FUNCTION  F2(Y) 

DOUBLE  PRECISION  Y,FACTRl,FACTR2,SOSQRD,S 
DOUBLE  PRECISION  SQRT2 
INTEGER  N 

COMMON  SOSQRD,S,N 

C  INITIALIZE  CONSTANTS 

SQRT2  =  DSQRT(2. ODO) 

C  COMPUTE  F2(Y,S,N) 

FACTR1  =  DERF(  Y/SQRT2  )**(N-2) 

FACTR2  =  DERF(  (Y+S)/SQRT2  ) 

F2  =  FACTRl*FACTR2*DEXP(-(Y**2)/2) 

RETURN 

END 


****************************************************************?v****>v** 
DOUBLE  PRECISION  FUNCTION  F3(Y) 

DOUBLE  PRECISION  Y.FACTRl .FACTR2 ,S0SQRD,S 
DOUBLE  PRECISION  SQRT2 
INTEGER  N 


COMMON  SOSQRD.S ,N 


C 

SQRT2  =  DSQRT(2.  ODO) 


INITIALIZE  CONSTANTS 


C  COMPUTE  F3(Y,S ,N) 

FACTR1  =  DERF(  (Y+S)/SQRT2  )**(N-2) 

FACTR2  =  DERF(  Y/SQRT2  ) 

F3  =  FACTRl*FACTR2*DEXP(-(Y+S)**2/2  ) 


RETURN 

END 
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